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A  USER'S  GUIDE  TO  NTRM,  A  SOIL-CROP  SIMULATION  MODEL  FOR 
NITROGEN,  TILLAGE,  AND  CROP -RESIDUE  MANAGEMENT 

M.J.  Shaffer  and  F.J.  Pierce1 


INTRODUCTION 

The  purpose  of  this  manual  is  to  provide  the  user  with  detailed 
operating  instructions  for  the  NTRM  model  (figs.  1  and  2).   This 
is  not  a  programmer's  manual,  and  you  should  not  expect  to  make 
model  changes  or  additions  based  on  the  information  given  here, 
Before  using  this  manual,  read  the  related  references  pertaining 
to  the  model's  concepts  and  theoretical  basis  (Shaffer  et  al . 
1983,  Shaffer  and  Larson  1987).   The  manual  includes  specific 
discussions  of  the  types  of  data  needed  to  run  the  model.   Use 
or  generation  of  all  input  and  output  media  such  as  printed 
output,  disk  files,  and  plots  are  discussed  in  detail  and  a 
general  output  file  summary  is  presented  in  appendix  A.   Use  of 
NTRM  is  described  for  each  major  part  of  the  model.   This 
includes  a  discussion  of  the  interactive  programs  that  are 
available  to  create  input  files  for  NTRM.   It  is  possible  to 
enter  the  required  input  data  and  then  receive  output  in  the 
form  of  screen  listings,  printouts,  or  graphics.   Sample  runs 
have  been  included  as  examples  of  how  to  use  the  model.   These 
are  not  exhaustive,  and  the  casual  user  may  want  to  consult  the 
detailed  NTRM  model  applications  manual  being  developed  to 
illustrate  the  many  uses  for  the  model. 

The  source  program  is  written  in  FORTRAN  77  computer  language 
for  Control  Data  Corporation  CYBER,  CRAY  Research  CRAY-1A,  HP 
A700,  and  IBM  PC/XT/AT  machines  or  true  compatibles.   The 
program  as  compiled  under  the  IBM  Professional  Fortran  (PROFORT) 
compiler  requires  about  512K  to  run  on  IBM  PC/XT /AT  machines. 

Units  in  this  manual  are  mixed  between  metric  and  English. 
Metric  units  are  used  except  where  experience  has  demonstrated  a 
greater  demand  for  the  English  equivalent.   Future  versions  of 
NTRM  will  contain  an  option  to  allow  the  user  to  select  the 
units  that  best  fit  local  requirements. 
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Figure  1 

Nitrogen-Tillage-Residue-Management 
(NTRM)  model. 
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Figure  2 

Processes  simulated  by  the  NTRM  model 

(with  saturated  zone). 


NTRM  MODEL  INSTALLATION,  SYSTEM 
CONFIGURATION,  AND  OPERATION 


The  NTRM  model  is  available  on  floppy  diskettes  which  contain  a 
copy  of  the  code  for  the  major  submodels  and  operating  system 
together  with  a  copy  of  sample  input  data  and  corresponding 
output.   These  programs  run  on  IBM  PC,  XT,  and  AT  computers  or 
true  compatibles  with  at  least  512K  RAM  (640K  recommended) ,  a 
hard  disk  drive,  floating  point  processor  8087  or  equivalent, 
and  DOS  2 . 1  or  newer.   The  FORTRAN  77  source  code  has  been 
included  along  with  the  object  files  for  each  program.   In  most 
cases,  the  linked  executable  files  have  also  been  included.   It 
should  be  possible  to  compile  and  link  the  source  code  on 
computer  systems  other  than  IBM  PC/XT/AT  provided  they  are 
compatible  with  FORTRAN  77. 

To  get  started  with  the  model  on  the  IBM  PC/XT/AT  or  true 
compatible,  we  suggest  that  you  set  up  a  separate  NTRM 
subdirectory  on  your  hard  disk  and  copy  all  the  files  from  the 
supplied  diskettes  into  that  directory.   You  should  also  modify 
your  CONFIG.SYS  file  to  allow  18  files  and  18  buffers.  The 
linked  NTRM  model  is  too  large  to  fit  on  a  single  360K  floppy 
diskette  and  must  be  stored  on  a  hard  disk  or  a  1.2  Mbyte  floppy 
diskette.   The  object  files  should  be  linked  using  IBM  Linker 
version  2.3  or  newer. 

Run  the  LKNTRMDE.BAT  linker  command  file  from  the  hard  disk 
after  the  object  files  have  been  placed  in  the  NTRM 
subdirectory.   Let  the  linker  create  the  default  files  by 
pressing  RETURN  when  prompted.   The  executable  file  will  be 
called  NTRM1DE.EXE.   In  a  similar  manner,  you  should  link  the 
soil  temperature  submodel  by  entering  LINK  NTRMST.OBJ  and 
pressing  RETURN.   The  executable  file  for  the  soil  temperature 
submodel  will  be  called  NTRMST.EXE.   Test  run  the  NTRM  main 
model  by  entering  NTRM1DE.EXE  and  pressing  RETURN.   Depending  on 
which  machine  you  have,  the  program  will  run  to  completion  in 
about  40  to  60  minutes.   Compare  the  contents  of  the  output 
files  just  generated  with  the  test  files  supplied  (same  names 
but  with  .TES  extension  added).   In  a  similar  manner,  test  the 
soil  temperature  submodel  using  NTRMST.EXE  and  compare  the 
output  created  on  files  SOTEMP  and  AIRTEMP  with  SOTEMP.TES  and 
AIRTEMP.TES.   Finally,  test  the  NTRM  data  file  setup  programs  by 
entering  RTEMP2.EXE  and  RNTRM2.EXE  followed  by  RETURN. 

The  general  configuration  of  the  NTRM  model  is  shown  in  figure 
3.   The  entire  model  can  be  operated  with  user- supplied  data  by 
first  using  the  appropriate  input  data  setup  routines  and  then 
running  the  appropriate  executable  files.   It  is  also  possible 
to  obtain  plots  of  data  generated  in  current  or  previous  model 
runs.   Note  that  model  input  is  file-oriented  and  that  certain 
input  files  must  be  present  in  correct  format  whenever  the  model 
is  run.   The  soil  temperature  submodel  is  a  separate  program 
with  its  own  required  inputs .   These  include  max-min  air 
temperatures  and  a  general  input  file  containing  control  data, 
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Figure    3 

NTRM  model  —  general  input  and 

output  configuration. 

model  parameter  values,  and  initial  conditions.   This  submodel 
is  normally  run  first,  and  its  output  becomes  input  to  the  NTRM 
main  interactive  model. 


This  central  soil-crop  model  requires  input  files  containing  (1) 
the  water  applications;  (2)  pan  evaporation,  wind,  and  solar 
radiation;  and  (3)  general  inputs  for  control  data,  initial 
conditions,  and  model  parameters. 

File  names  for  running  the  central  model  must  be  established  in 
a  file  called  NTRFILE.   Names  for  all  required  and  optional 
input  files  and  optional  output  files  are  specified  in  this 
file.   File  names  are  limited  to  seven  characters  or  less. 
Output  will  be  created  in  files  with  the  names  specified.   If  a 
particular  output  file  is  not  desired,  the  word  NONE  should  be 
entered  and  the  corresponding  file  will  not  be  opened  or 
written.   The  following  input  and  output  file  names  must  be 
specified  in  NTRFILE  in  the  order  shown  below  (one  name  per 
line) . 


Contents  of  NTRFILE: 

Line  No.  Description  Required  (Yes  or  No) 

1 .  Common  input  file  Yes 

2 .  General  output  file  No 

3.  Soil  temperature  input  file  No 

4.  Air  temperature  input  file  No 

5 .  Precipitation  and  irrigation 

input  file  Yes 

6.  Pan  evaporation,  solar,  and 

wind  input  file  Yes 

7.  Deep  percolation  output  file  No 

8 .  Top  growth  output  file  No 

9 .  Root  dry  matter  output  file  No 

10.  Daily  roots  summary  No 

11.  Unsaturated  flow  output  No 

12 .  Nitrogen  output  file  No 

13.  3-D  plot  output  file  No 

14 .  Surface  data  output  file  No 

15.  N15  output  file  No 

16.  Evapotranspiration  output  file  No 

17 .  Special  graphics  control  file  No 
18  .  Unsaturated  flow  data  (may  be 

input  or  output)  No 

The  aquifer  submodels  use  output  from  the  soil -crop  model  and 
require  additional  general  inputs .   These  submodels  have  been 
described  previously  by  Ribbens  (1976)  and  by  Shaffer  and  Gupta 
(1981) . 

Plot  files  generated  by  the  various  submodels  are  in  a  form 
suitable  for  listing  or  direct  input  to  most  plotting  packages. 
A  special  interactive  graphics  system  is  available  from  the 
authors  which  allows  direct  plotting  of  NTRM  plot  files  on 
TEKTRONIX  and  other  graphical  equipment.   This  system  is  a 
separate,  stand-alone  package  that  supplements  the  NTRM  software 
described  in  this  manual.   In  general,  output  from  NTRM  is 
suitable  for  input  to  a  line  printer;  exceptions  are  certain 
optional  binary  files  such  as  one  version  of  the  deep 
percolation  file . 

The  NTRM  interactive  input  file  setup  programs  simplify  the  use 
of  the  model  on  microcomputers  such  as  IBM  PC/XT/AT .   These 
programs  were  developed  to  simplify  application  of  the  model. 
Batch  and  card  deck  operation  are  also  still  possible  via  a 
mainframe  or  minicomputer,  but  the  user  is  responsible  for 
installation  (including  compilation)  and  maintenance  of  NTRM  on 
these  systems.   Each  file  setup  program  is  entered  with  a  single 
system  command.   A  series  of  menus  and  questions  are  then 
displayed  to  guide  you  through  creation  of  the  required  input 
files . 


SAMPLE  PROBLEM 


One  of  the  best  ways  to  learn  how  to  operate  a  complex  model 
such  as  NTRM  is  to  follow  through  an  example.   For  this 
exercise,  we  will  assume  that  you  wish  to  make  runs  using  your 
IBM  PC/XT/AT  or  true  compatible  machine.   Your  microcomputer 
should  be  equipped  with  a  printer  to  record  the  data  entered  and 
the  results.   Graphics  capabilities  in  the  form  of  a  display  . 
screen  and/or  plotter  are  highly  desirable,  particularly  if  the 
equipment  is  compatible  with  TEKTRONIX  code. 

For  this  example,  we  will  assume  that  you  wish  to  make  a  short 
test  run  with  the  NTRM  model  using  data  contained  in  a  field 
notebook.   The  most  common  starting  place  for  NTRM  runs  is  to 
create  two  files  containing  (1)  the  daily  max-min  air 
temperatures  and  (2)  the  general  input  data  for  the  soil 
temperature  submodel.   This  can  be  done  by  running  the  NTRM  file 
set-up  program  for  soil  temperature  called  RTEMP2.EXE.   In 
addition,  data  could  be  uploaded  from  a  remote  microcomputer 
disk  or  tape  unit.   The  data  formats  for  these  two  files  are 
described  in  "Common  Input  File  for  Soil  Temperature  Submodel" 
and  in  "Max-Min  Air  Temperature  Input  File  for  Soil  Temperature 
Submodel".   Once  the  files  are  saved  on  the  disk,  the  NTRM  soil 
temperature  submodel  NTRMST.EXE  can  be  used  to  make  the  actual 
soil  temperature  simulation  runs. 

Please  note  that  the  microcomputer  version  of  the  soil 
temperature  submodel  requires  that  the  daily  max-min  air 
temperature  file  be  named  AIR,  and  the  general  temperature  input 
data  file  be  named  INPUT1 .   Once  the  AIR  and  INPUT1  files  are 
prepared,  the  soil  temperature  simulation  is  run  by  typing 
NTRMST.EXE  and  pressing  the  RETURN  key. 

Please  refer  to  appendix  B  for  the  results  to  expect  when  the 
NTRM  operating  system  is  used  to  set  up  the  necessary  soil 
temperature  files  and  make  a  soil  temperature  simulation  run. 
In  this  example,  the  max-min  air  temperatures  are  taken  from  the 
notebook  and  entered  directly  via  the  keyboard  of  the 
microcomputer.   The  general  input  file  is  also  created  using  the 
same  program.   You  are  required  to  supply  only  initial  soil 
temperatures,  Julian  starting  date,  length  of  the  run,  and 
management  code.   Other  parameters  are  set  to  default  values  by 
the  operating  system.   These  can  be  modified  by  using  one  of  the 
text  editors  on  your  microcomputer,  such  as  EDLIN. 

This  particular  simulation  run  created  daily  soil  temperature 
and  hourly  air  temperature  files  named  SOTEMP  and  AIRTEMP , 
respectively.   These  are  suitable  for  direct  input  into  the  rest 
of  the  NTRM  model.   The  output  file  named  OUTPUT  contains 
general  information  in  tabular  form  for  listing  on  a  line 
printer. 


The  next  step  in  making  the  overall  NTRM  run  involves  creating  a 
second  general  input  file  (the  common  input  file),  a  water 
application  file,  and  a  pan  evaporation  file  with  the  help  of 
the  NTRM  file  setup  program  RNTRM2.EXE.   Refer  again  to  appendix 
B  for  the  results  obtained  from  a  sample  computer  session.   Note 
that  each  of  the  required  files  can  be  created  and  stored  on  the 
disk  using  the  NTRM  operating  system.   In  this  case,  the  general 
input  file  will  always  contain  some  variables  that  are  set  to 
default  values.   The  MS-DOS  text  editor  can  be  used  to  modify 
these  if  necessary.   It  is  also  possible  to  create  these  files 
off  line  with  a  microcomputer  or  similar  device  and  then  upload 
them  to  your  machine.   In  addition,  the  files  can  be  entirely 
created  on  your  microcomputer  by  using  EDLIN  or  some  other  text 
editor  or  word  processor  program.   This  often  involves 
modification  of  an  existing  file(s). 

Once  these  files  are  available  on  your  microcomputer  disk,  the 
NTRM  program  named  NTRM1DE.EXE  can  be  run  to  simulate  the 
remainder  of  the  soil-plant  processes.   However,  first  be 
certain  that  all  the  file  names  have  been  correctly  specified  in 
the  control  file  NTRFILE.   Once  this  has  been  done,  the 
simulation  is  initiated  by  typing  NTRM1DE.EXE  and  pressing 
RETURN.   Typical  execution  times  on  an  IBM  AT  computer  for  a 
crop-growing  season  range  from  20  to  50  minutes  depending  on  the 
length  of  the  run. 

It  is  possible  to  save  and  preview  plot  files  before  the  actual 
graphics  output  files  are  created.   This  often  saves  time  and 
money  when  a  new  set  of  input  data  is  being  tried.   This  is  done 
by  specifying  in  NTRFILE  which  plot  files  to  create  and  then 
bypassing  (temporarily)  the  plots  themselves.   The  plot  files 
can  be  listed  at  a  terminal  or  line  printer  for  review.   These 
files  can  then  be  converted  to  graphics  output  files  by  calling 
an  appropriate  graphics  package.   It  is  the  graphics  output 
files  that  are  suitable  for  direct  input  into  TEKTRONIX- 
compatible  display  tubes  or  plotters  or  into  line  printers  (for 
plots) . 


INTERACTIVE  GENERATION  OF  INPUT  FILES 


This  section  describes  the  portion  of  the  NTRM  model  that 
assists  you  in  interactive  input  of  soil  chemical  and  physical 
properties,  crop  growth  data,  management  data  (tillage  and 
fertilizer),  and  climate  data  necessary  to  run  the  model.   These 
data  inputs  are  discussed  in  the  sections  of  the  manual  on  the 
operation  of  each  submodel.   When  you  operate  the  interactive 
system,  please  remember  that  some  of  the  parameter 
specifications  are  set  to  default  values  and  are  therefore  not 
immediately  available  for  user  input.   The  number  of  parameters 
in  this  category  varies  with  the  input  option  selected.   You  can 
always  modify  any  of  the  parameter  values  later  by  using  the 
EDLIN  or  similar  text  editor. 

When  you  run  either  RTEMP2.EXE  or  RNTRM2.EXE,  the  operating 
system  will  ask  you  to  enter  names  of  files  you  wish  to  create. 
With  the  soil  temperature  file  program  RTEMP2.EXE,  the  system 
responds 

ENTER  CONTROL  CODE: 

1  =  CREATE  MAX-MIN  AIR  TEMPERATURE  FILE 

2  =  CREATE  GENERAL  SOIL  TEMPERATURE  INPUT  FILE 

3  =  CREATE  BOTH  FILES 

and  with  the  central  NTRM  model  file  program  RNTRM2.EXE,  the 
system  prints 

ENTER  FILES  TO  CREATE: 

1  =  COMMON  INPUT  FILE 

2  =  WATER  APPLICATION  FILE 

4  =  PAN  EVAPORATION,  WIND,  AND  SOLAR  RADIATION  FILE 
7  =  ALL  OF  THE  ABOVE 

NOTE:   TO  GET  COMBINATIONS  OF  THE  ABOVE,  SUM  THE 
NUMBERS  INVOLVED 

Any  combination  of  these  files  can  be  created.   For  each 
selection,  you  will  be  asked  to  enter  a  filename.   You  are  also 
offered  the  option  of  stopping  the  program  after  each  file 
creation  or  of  continuing  to  the  next  portion  of  the  operating 
system. 

Remember,  the  operating  system  can  be  stopped  at  any  time  by 
typing  CTRL  C  in  response  to  any  input  queue.   This,  however, 
may  result  in  loss  of  all  data  input  in  that  procedure.   A 
complete  set  of  sample  runs  is  presented  in  appendix  B.   These 
demonstrate  the  application  of  the  model  to  an  actual  problem. 


Creating  Max-Min 
Air  Temperature 
File 


This  subprogram  asks  you  to  enter  the  Julian  starting  day  and  a 
series  of  daily  maximum  and  minimum  air  temperatures .   Enter  the 
temperature  values  free  format  with  the  maximum  coming  before 
the  minimum  in  each  case .   Free  format  means  that  one  or  more 
values  separated  by  blanks  or  commas  can  be  entered  anywhere  on 
the  same  line.   Stop  entry  of  air  temperature  data  by  pressing 
RETURN  in  response  to  the  prompt  for  a  new  line  of  data. 


Creating  General 
Soil  Temperature 
Input  File 


This  routine  allows  you  to  rapidly  create  a  general  file  for  use 
in  making  soil  temperature  simulation  runs.   You  supply  values 
for  initial  soil  temperatures,  key  control  parameters,  and 
management  options.   Remaining  parameter  values  are  provided  by 
default.   These  can  be  altered,  if  necessary,  by  using  one  of 
the  text  editors.   Consult  the  appropriate  sections  of  this 
manual  for  formatting  details . 


Creating  Water 
Application  File 


The  flow  of  data  input  for  this  file  is  shown  in  figure  4.   The 
program  to  create  the  water  application  file  has  two  parts. 
First,  you  are  asked  to  input  the  number  of  water  applications 
(up  to  99)  and  some  coded  information. 


THIS  PROGRAM  GUIDES  THE  USER  IN  THE 
INPUT  OF  WATER  APPLICATIONS  DATA  FOR  NTRM 

ENTER  A  NUMERIC  IDENTIFIER  CODE 

ENTER  YEAR  OF  THE  RUN,  EX.,  1981 

ENTER  NUMBER  OF  WATER  APPLICATIONS 
(MAXIMUM  NUMBER  IS  99) 

ENTER  ALPHANUMERIC  DESCRIPTION  TO  IDENTIFY 
THE  WATER  APPLICATIONS 

Next  you  are  prompted  to  input  the  data  and  type  of  water 
application  (R  =  rain;  I  =  irrigation)  and  the  amount  (vol/unit 
area)  for  each  water  application. 

ENTER  MONTH,  DATE  OF  MONTH, 

AND  AN  ALPHANUMERICAL  IDENTIFIER  CODE 

FOR  THE  TIME  OF  WATER  APPLICATION  1  AND  TYPE 

(EX.,  0121R=JANUARY  21,  RAIN) 

ENTER  AMOUNT  OF  APPLICATION  1  AS  A 
VOLUME  PER  UNIT  AREA  (EX. ,  INCHES) 

On  completion  of  input  for  the  last  water  application,  program 
control  returns  to  the  operating  system,  which  selects  the  next 
procedure  to  run. 
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f  START         J 

/enter  file/ 

/    NAME   / 


ENTER  #       | 
APPLICATIONS  AND  / 
CODING  INFORMATION/ 


FOR  1=1  TO 
#  APPLICATIONS 


f 


CgotonextN 
procedure/ 


NO 


Figure  4 

Flow  of  data  input  for  the  creation  of 

water  application  file. 
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Creating  Pan 
Evaporation,  Wind, 
and  Solar  Radiation 
File 


The  flow  of  the  data  input  for  this  file  is  shown  in  figure  5. 
Two  options  are  presently  available  to  enter  pan  evaporation 
data- -semimonthly  or  daily  (+  wind  travel  and  solar  radiation) 

THIS  PROGRAM  GUIDES  THE  USER  IN  THE  INPUT 
OF  PAN  EVAPORATION  DATA  EITHER  ON  A  SEMI-MONTHLY  BASIS 
OR  A  DAILY  BASIS.   IN  THE  LATTER  CASE,  DAILY  WIND 
TRAVEL  (KM/DAY)  AND  SOLAR  RADIATION  (LANGLEYS)  ARE 
ALSO  REQUIRED. 


Figure  5 

Flow  of  data  input  forthe  creation  of 
pan  evaporation,  wind  run,  and  solar 
radiation  file. 
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CHOOSE  ONE  OF  THE  FOLLOWING  OPTIONS  FOR  PAN  EVAPORATION 

0.  SEMI-MONTHLY  PAN  EVAPORATION 

1.  DAILY  PAN  EVAPORATION  (+  WIND  +  SOLAR) 

Semimonthly 

For  this  option  you  must  input  a  numeric  code,  year  of  data,  and 
24  semimonthly  pan  evaporation  data  (inches). 

ENTER  A  NUMERIC  CODE  (MAX  5  CHARACTERS)  THAT  CAN  BE  USED 
TO  IDENTIFY  THE  CROP,  CLIMATE,  IRRIGATION  METHOD,  ETC. 
(DISPLAYED  ON  OUTPUT) 

ENTER  YEAR  OF  THE  RUN,  EX. ,  1981 

FOR  EACH  OF  24  PERIODS  FOR  THE  YEAR,  ENTER  THE  SEMI-MONTHLY 
PAN  EVAPORATION  IN  INCHES 

ENTER  PAN  DATA  FOR  PERIOD  1 

Daily 

For  this  option,  you  must  enter  the  number  of  days  of  the 
simulation  run  and  the  pan  evaporation  (inches)  and  wind  travel 
(km/day)  and  solar  radiation  (langleys/day)  for  each  day  of  the 
run. 

THE  PROGRAM  WILL  PROMPT  YOU  IN  THE  INPUT 

OF  DAILY  PAN  EVAPORATION,  WIND  TRAVEL,  AND 

SOLAR  RADIATION.   THE  PAN  EVAPORATION  IS  IN  UNITS 

OF  INCHES,  WIND  IS  IN  UNITS  OF  KM/DAY, 

AND  SOLAR  RADIATION  IS  IN  UNITS  OF  LANGLEYS/DAY. 

(EX. ,  .217,  193,  464) 

ENTER  THE  NUMBER  OF  DAYS  OF  THE  SIMULATION  RUN 
ON  THE  NTRM  MODEL 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY  1 

When  all  values  have  been  entered,  program  control  will  return 
to  the  operating  system,  which  selects  the  next  procedure  to 
use . 
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Creating  Common 
Input  File  for 
Central  Model 


The  flow  of  data  input  for  this  file  is  shown  in  figure  6 
system  first  requests  the  selection  of  an  input  option. 


The 


SELECT  SUBMODEL 
CONTROL  AND 
PRINT  OPTIONS  (2,25) 


TITLE  AND  TIME 
FRAMEWORK  INPUT  (1) 


SOIL  HORIZONS 
AND  DEPTH  (A) 


I 


#  SEGMENTS  FOR 
EACH  HORIZON 


(5)1 


CHEMICAL  DATA  FOR 
EACH  HORIZON  (7&13) 


Figure  6 

Flow  of  data  input  for  the  creation  of 
common  input  file  for  central  NTRM  model. 
Numbers  in  parentheses  correspond  to 
data  group. 


1/. 


PHYSICAL  DATA  FOR 
EACH  HORIZON  (26) 


ENTER  EVAPORATION 
MODEL  DATA  (24) 


f 


rENTER  MOISTURE 
FLOW  DATA  FOR 
EACH  SOIL 
SEGMENT 


INPUT 
TEMPERATURE 
HORIZONS 
AND  WEEKLY 
TEMPERATURE 


SET  RAIN  AND/OR 
IRRIGATION  WATER 
ANALYSIS 
DEFAULT  =0.0  (9,10) 


Figure  6 — Continued 


L5 


CROP  GROWTH 
DATA  (23) 


TILLAGE  DATES,  DEPTHS 
AND  TYPE  (17) 


INPUT  #  FERTILIZER 
(INORGANIC) 
APPLICATIONS 
(15,19) 


APPLICATION  DATES 
DEPTH  OF  APPLICATION 
CONTENTS  (LBS / AC ) 


INPUT  #  ORGANIC 
-V  APPLICATIONS 
(16,20) 


APPLICATION  DATES 
DEPTH  OF  APPLICATION 
CONTENTS   (LBS/ AC) 


RETURN  TO 

OPERATING 

SYSTEM 


Figure  6 — Continued 
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GENERATION  SYSTEM  FOR  COMMON  INPUT  FILE 

THE  FOLLOWING  ARE  THE  INPUT  OPTIONS  PRESENTLY  AVAILABLE 

IN  THIS  OPERATING  PROCEDURE 

1.  DEFAULT  OPTION  1  --  MINIMAL  USER  INPUT  REQUIRED 

2.  DEFAULT  OPTION  2  --  ALLOWS  SELECTION  OF  SUBMODELS 

UTILIZED,  MORE  USER  INPUT  MAY 
BE  REQUIRED 

ENTER  OPTION  SELECTION 

Option  1  is  the  short  option,  which  minimizes  the  inputs 
required  to  run  the  model.   This  option  is  illustrated  in  the 
sample  output  presented  in  Appendix  B.   Option  2  expands  the 
number  of  parameters  you  can  set,  including  submodel,  control, 
and  print  option  selection.   If  option  2  is  selected,  the  status 
of  the  submodel  and  control  options  are  listed  at  the  terminal 
and  the  user  is  asked  to  make  the  appropriate  changes.   When  all 
the  changes  are  complete,  the  new  status  of  these  options  is 
again  listed  at  the  terminal.   The  status  of  the  print  options 
is  then  listed  at  the  terminal,  and  you  are  asked  to  make  the 
appropriate  changes  in  these  options. 


THE  FOLLOWING  IS  A  SUMMARY  OF  THE 
OPTIONAL  SUBMODELS  AND  CONTROL  SWITCHES 

OPTIONAL  SUBMODEL  OR  SWITCH  DEFAULT 

NITROGEN  TRANSFORMATIONS  0 

A.  0  IF  CONSIDERED 

B.  1  IF  BYPASSED 

SALT  REACTIONS  1 

A.  0  IF  CONSIDERED 

B.  1  IF  BYPASSED 

UNSATURATED  FLOW  MODEL  1 

A.  0  IF  READ  FROM  TAPE  1  OR  TERMINAL  INPUT 

B.  1  IF  UNSATURATED  FLOW  SUBMODEL  UTILIZED 

SOIL  TEMPERATURE  MODEL  1 

A.  0  IF  WEEKLY  SOIL  TEMPERATURES  INPUT 

B.  1  IF  SOIL  TEMPERATURE  SUBMODEL  UTILIZED 

EVAPORATION  MODEL  1 

A.  0  IF  POTENTIAL  EVAPORATION  IS  SET  TO  ZERO 

B.  1  IF  POTENTIAL  EVAPORATION  IS  ESTIMATED 
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6.  WATER  FLOW  DATA  0 

A.  0  IF  TO  BE  READ  FROM  TAPE  1 

B.  1  IF  TO  BE  INPUT  FROM  TERMINAL 

7.  DEMAND  FOR  CROP  GROWTH  0 

A.  0  IF  CROP  GROWTH  INTERACTIVE  WITH 
NITROGEN  AVAILABILITY 

B.  1  IF  CROP  GROWTH  IS  IDEAL  WITH  RESPECT 
TO  NITROGEN 

8.  ROOT  GROWTH  CONTROL  SWITCH  0 

A.  0  IF  ROOT  GROWTH  IS  INTERACTIVE  WITH 
SOIL  CONDITIONS 

B.  1  IF  ROOT  GROWTH  IS  IDEAL  WITH  RESPECT 
TO  SOIL  CONDITIONS 

DO  YOU  WISH  TO  CHANGE  ANY  OF  THE  ABOVE  OPTIONS? 

ENTER  1  FOR  YES,  0  FOR  NO 

?  1 

ENTER  OPTION  NUMBER  YOU  WISH  TO  CHANGE 
(IF  CHANGES  ARE  COMPLETE,  ENTER  0) 

?  2 

ENTER  NEW  VALUE  (1  OR  0)  FOR  OPTION  2 

?  1 

ENTER  OPTION  NUMBER  YOU  WISH  TO  CHANGE 
(IF  CHANGES  ARE  COMPLETE,  ENTER  0) 

?  0 

Some  of  the  submodel  options  are  not  available  using  the  above- 
mentioned  switches  but  rather  are  selected  by  the  presence  of 
additional  data  that  are  added  to  the  common  input  file.   These 
data  must  be  added  to  an  existing  common  input  file  for  the 
central  NTRM  model. 

THE  FOLLOWING  IS  THE  CURRENT  STATUS  OF  THE  OUTPUT  OPTIONS 

1.  PRINT  INTERVAL  (DAYS)  FOR  PREDICTED  SOIL  CHEMISTRY      2 

2.  PRINT  INTERVAL  (DAYS)  FOR  SUMMARY  OF  NITROGEN 

BALANCE  FOR  SYSTEM  10 

3.  TIME  INTERVAL  WITHIN  A  DAY  ON  WHICH  1.  AND  2.  ABOVE 

ARE  TO  BE  ACTIVATED  10 


4.  WRITE  INTERVAL  (DAYS)  FOR  OUTPUT  TO  TAPE  2  2 

5.  PRINTING  OF  BASIC  CONTROL  CARD  DATA  1 

6.  PRINTING  OF  INPUT  DATA  1 

7.  WRITE  N15  DATA  TO  OUTPUT  TAPE  15  0 

8.  PRINTING  OF  INITIAL  MOISTURE  INPUT  DATA  1 

9.  PRINTING  OF  COMPUTED  MOISTURE  CONTENTS  ON  A 

DAILY  BASIS  1 

DO  YOU  WISH  TO  CHANGE  ANY  OF  THE  PRINT  OPTIONS  LISTED 
ENTER  1  FOR  YES,  0  FOR  NO 
?  0 

Title  and  Time  Framework 

You  are  prompted  to  enter  output  title,  beginning  and  ending 
date  (Julian)  of  the  run,  and  the  date  of  crop  emergence.   If 
you  have  selected  the  plotting  option,  the  title  should  end  with 
the  character  "$." 

INPUT  TITLE  DESIRED  ON  PRINTED  OUTPUT 
(1  TO  80  CHARACTERS) 

THE  FOLLOWING  INPUT  SETS  TIME  FRAMEWORK 
FOR  THE  SIMULATION  RUN 
INPUT  JULIAN  DATES  FOR  THE  FOLLOWING 
STARTING  DATE  OF  THE  RUN  = 

TERMINATION  DATE  OF  THE  RUN  = 

ENTER  JULIAN  DATE  FOR  CROP  EMERGENCE 


Soil  Horizons ,  Depths ,  and  Segments 

You  are  asked  to  input  the  number  of  soil  horizons  (up  to  10) 
and  the  depth  (cm)  to  the  lower  horizon  boundary  for  each 
horizon.   Maximum  depth  of  the  soil  profile  is  not  allowed  to 
exceed  150  cm. 

THE  FOLLOWING  INPUT  SETS  SOIL  DATA 
INPUT  NUMBER  OF  SOIL  HORIZONS  (1  TO  10) 


INPUT  DEPTH  (CM)  FROM  SURFACE  TO  BOTTOM  OF  HORIZON  1 
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This  last  prompt  is  repeated  for  each  soil  horizon. 


c 
: 
: 


Next,  you  are  asked  to  enter  the  number  of  segments  per  horizon. 

We  recommend  that  the  total  number  of  segments  for  the  entire 

soil  profile  be  at  least  10  to  preclude  error  in  the  dispersion 

simulation  (a  message  to  this  effect  will  be  printed  at  the 

terminal  if  fewer  than  10  segments  are  included) .  m 

ENTER  NUMBER  OF  SEGMENTS  FOR  HORIZON  1 


ENTER  NUMBER  OF  NODAL  SUBDIVISIONS  IN 
EACH  SEGMENT  (FOR  HORIZON  1) 


This  will  be  repeated  for  each  horizon. 

Soil  Chemical  Data  (Data  Groups  7  and  13) 

Both  options  1  and  2  require  input  of  soil  chemical  information. 

The  required  input  for  option  1  for  each  chemistry  horizon  is         <- 

marked  with  an  *  in  the  example  output  (see  below) .   The  amount 

of  input  required  for  option  2  depends  on  the  submodel  section 

for  the  nitrogen  transformation  and  salt  reactions .  — 

1.  If  both  the  nitrogen  transformations  and  the  salt  „, 
reactions  are  not  considered  (status  =  0) ,  then  only  the 

bulk  density  for  each  horizon  is  required  input. 

m 

2.  If  the  nitrogen  transformations  are  considered  and  the 

salt  reactions  are  not  considered,  then  the  inputs  - 

required  are  the  same  as  those  for  option  1  (the  * 

inputs) .  f 

3.  If  the  nitrogen  transformations  are  not  considered  and  the 

salt  reactions  are,  then  those  inputs  marked  with  a  +  are        _ 
required. 

m 

4.  If  both  the  nitrogen  transformations  and  the  salt 
reactions  are  considered,  then  all  the  inputs  listed  in 

the  example  output  are  required.  L 

*  INPUT  THE  FOLLOWING  SOIL  CHEMICAL  ANALYSIS  FOR  CHEMICAL  - 
HORIZON  1 

*  ENTER  NH4+  (SOLUBLE  FORM)  (MEQ/L) 
(WARNING! !  NEVER  USE  TOTAL  NH4+  HERE) 
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: 


?  .027 

*  ENTER  N03"  (MEQ/L) 
?  .009 

*  ENTER  UREA  (MEQ/L) 
?  0 

ENTER  CA++  (MEQ/L) 

?  .2 

ENTER  NA+  (MEQ/L) 

?  .15 

ENTER  MG++  (MEQ/L) 

?  .003 

ENTER  HC03'  (MEQ/L) 

?  0 

ENTER  CL"  (MEQ/L) 

?  .02 

ENTER  C03=  (MEQ/L) 

?  .03 

ENTER  S04=  (MEQ/L) 

?  .03 

ENTER  EXCHANGE  CAPACITY  (MEQ/100G  SOIL) 

?  15 

ENTER  GYPSUM  (MEQ/100G  SOIL) 

?  0 

IS  LIME  PRESENT  IN  SOIL?   0.0  IF  NO,  1.0  IF  YES 

?  0.0 
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*  ENTER  BULK  DENSITY  (G/CM3  OF  SOIL) 
?  1.35 

*  ENTER  GM  WATER/GM  SOIL  IN  SOIL  EXTRACT 

?  5 

SINCE  SALT  REACTIONS  ARE  TO  BE  CONSIDERED,  ENTER  C02  PARTIAL 
PRESSURE  (ATM)  FOR  HORIZON  1 

?  .008 

*  ENTER  MICROBIAL  BIOMASS  (UG/G  SOIL)  AS  CARBON 
?  50 

*  ENTER  NUTRIENT  HUMUS  (UG/G  SOIL)  AS  CARBON 
?  50 

*  ENTER  MAXIMUM  SURFACE  INFILTRATION  RATE  (CM/DAY) 
?  50.0 

Soil  Temperature  (optional) 

If  option  1  is  selected,  you  must  have  the  soil  temperature  data 
on  file  in  your  catalog  and  this  section  of  input  will  be 
bypassed.   If  option  2  is  selected  and  you  set  the  soil 
temperature  parameter  equal  to  zero,  then  you  must  enter  the 
soil  temperature  horizons  and  weekly  temperatures  as  follows : 

THE  FOLLOWING  SETS  UP  TEMPERATURE  HORIZONS  AND  DATA 
ENTER  NUMBER  OF  TEMPERATURE  HORIZONS  (1  TO  4) 

?  1 

ENTER  DEPTH  (CM)  FROM  SURFACE  TO  BOTTOM  OF 
TEMPERATURE  HORIZON  1 

?  150 

ENTER  NUMBER  OF  WEEKLY  TEMPERATURES  FOR  THIS  RUN 

?  2 

ENTER  MEAN  SOIL  TEMPERATURE  (C)  FOR  WEEK  1 
FOR  TEMPERATURE  HORIZON  1 

?  15 
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ENTER  MEAN  SOIL  TEMPERATURE  (C)  FOR  WEEK  2 
FOR  TEMPERATURE  HORIZON  1 

?  25 

If  you  have  the  soil  temperature  on  file,  then  set  the  soil 
temperature  status  as  listed  in  the  status  section  equal  to  1 
and  this  section  of  input  will  be  bypassed.   There  is  an  option 
in  the  operating  system  to  run  the  soil  temperature  model. 

Waterflow  Data  (optional) 

For  option  1,  the  unsaturated  flow  model  is  selected  by  default. 
Therefore  this  section  is  bypassed.   For  option  2,  this  section 
is  not  necessary  if  the  unsaturated  flow  model  is  selected  or  if 
you  have  waterflow  data  on  tape.   Note  that  the  amount  of  input 
required  for  this  option  can  be  large. 

THE  FOLLOWING  GUIDES  THE  USER  IN  THE  INPUT 
OF  WATER  FLOW  DATA  WHEN  THE  UNSATURATED 
FLOW  MODEL  IS  NOT  USED  AND  THE  FLOW  DATA 
IS  NOT  ALREADY  ON  A  DATA  TAPE. 


ENTER  WATER  FLOW  DATA  FOR  SOIL  SEGMENTS 

ENTER  DEPTH  (CM)  OF  WATER  ON  THE  SOIL  SURFACE 

?  1 

ENTER  WATER  FLUX  ACROSS  SOIL  SURFACE  (CM/DAY) 

?  3 

ENTER  AVERAGE  VOLUME  (CM  IN  UNIT  CROSS -SECTION) 
OF  WATER  IN  SEGMENT  1 

?  5 

ENTER  WATER  FLUX  (CM/DAY)  ACROSS  LOWER  BOUNDARY 
OF  SOIL  SEGMENT  1 

?  3 

ENTER  AVERAGE  WATER  TENSION  (CM)  OVER  TIME 

STEP  FOR  SOIL  SEGMENT  1  (TENSIONS  SHOWN  AS  NEGATIVE) 

?  -300 
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ENTER  PLANT  WATER  UPTAKE  (CM/DAY)  FOR  SOIL  SEGMENT  1 

?  1.5  L 

ENTER  AVERAGE  VOLUME  (CM  IN  UNIT  CROSS  -  SECTION)  OF  - 

WATER  IN  SEGMENT  2 

?  5 

This  is  repeated  for  each  segment  entered  previously.  — 

Rain  and  Irrigation  Water  Analysis 

This  is  always  required.   You  have  the  option  of  setting  either       — 
analysis  to  a  default  value  of  0.0. 

THE  FOLLOWING  WILL  GUIDE  IN  THE  INPUT  OF  RAIN  AND  IRRIGATION 
WATER  ANALYSIS 

DO  YOU  WISH  TO  SET  RAIN  WATER  ANALYSIS 

TO  DEFAULT  VALUE  OF  0.0.   ENTER  1  FOR  YES,  - 

0  FOR  NO 

?  1 

IS  IRRIGATION  WATER  TO  BE  APPLIED  - 

ENTER  1  FOR  YES,  0  FOR  NO 

?  1  L 

DO  YOU  WISH  TO  SET  ANALYSIS  TO  DEFAULT  VALUE 

OF  0.0.   ENTER  1  FOR  YES,  0  FOR  NO  - 

?  0 

If  you  choose  to  input  water  analysis,  the  following  input 

requests  will  appear  at  the  terminal.  L 

ENTER  THE  FOLLOWING  IRRIGATION  WATER  ANALYSIS 
(UNITS  ARE  MEQ/L) 

ENTER  NH4+ 

?  .027 

ENTER  N03" 
?  .009 
ENTER  Ca++ 


2  4 


?  2.0 


ENTER 

Na++ 

?  .5 

ENTER 

Mg++ 

?  1.0 

ENTER 

HC03 

?  2.0 

ENTER 

Cl" 

?  1.0 

ENTER 

C03= 

?  0.0 

ENTER 

S04= 

?  0.5 

Crop  Type  and  Coefficients 

This  information  is  always  required. 

ENTER  THE  FOLLOWING  CROP  GROWTH  DATA 

ENTER  CROP  MATURITY  CLASS  FOR  CORN  (IN  DAYS) 
EX.  85  =  85  DAY  CORN,  110  =  110  DAY  CORN,  ETC. 

?  85 

ENTER  PLANT  POPULATION  (PLANTS/HA) 

?  56198 

ENTER  UP  TO  10  CHARACTER  NAME  OF  CROP  VARIETY 

?  CORN  TEST 
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Tillage  Data  Input 

The  following  is  required  for  both  options: 

THIS  PORTION  OF  THE  PROGRAM  GUIDES  THE  USER  IN  ENTERING  TILLAGE 
INFORMATION 

ENTER  NUMBER  OF  TILLAGE  EVENTS  (MAX=6) 

ENTER  JULIAN  DATE  OF  TILLAGE  EVENT  1 

ENTER  DEPTH  (CM)  OF  TILLAGE  FOR  EVENT  1 

ENTER  TILLAGE  CODE.   1  =  MOLDBOARD ,  2  =  CHISEL  PLOW,  3  =  DISK 

Fertilizer  Applications 

You  are  asked  to  enter  the  number  of  fertilizer  applications. 
If  there  were  no  applications,  enter  0  and  the  remainder  of  the 
input  is  bypassed.   If  fertilizer  was  applied,  you  are  asked  for 
the  Julian  application  date,  the  depth  of  application,  and  the 
content  of  the  fertilizer  for  each  fertilizer  application. 

THE  FOLLOWING  GUIDES  IN  FERTILIZER  DATA  INPUT 
INPUT  THE  NUMBER  OF  FERTILIZER  APPLICATIONS 

?  1 

ENTER  APPLICATION  DAY  NUMBER  1  (JULIAN  DATE) 

?  121 

THE  FOLLOWING  WILL  SET  UP  FERTILIZER  APPLICATIONS 
FOR  EACH  FERTILIZER  APPLICATION  DATE 

ENTER  DEPTH  (CM)  OF  A  UNIFORM  FERTILIZER  APPLICATION 
FOR  APPLICATION  DATE  NUMBER  1 

(IF  =  0.0  A  SURFACE  APPLICATION  IS  INDICATED) 

DEPTH  = 

?  4 

ENTER  LBS/ACRE  OF  NH4+  ADDED  IN  APPLICATION  1 
TYPICAL  RANGE  0-400  LBS/ACRE 

?  25 
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ENTER  LBS/ACRE  OF  N03-  ADDED  IN  APPLICATION  1  TYPICAL  RANGE 
0-1600  LBS/ACRE 

?  40 


ENTER  LBS/ACRE  OF  UREA  ADDED  IN  APPLICATION  1  TYPICAL  RANGE 
0-800  LBS/ACRE 

?  0 

ENTER  LBS/ACRE  OF  CA++  ADDED  IN  APPLICATION  1  RANGE  0-20,000 
LBS/ACRE 

?  0 

ENTER  LBS/ACRE  OF  S04=  ADDED  IN  APPLICATION  1  RANGE  0-20,000 
LBS /ACRE 

?  50 

ENTER  LBS/ACRE  OF  C03=  ADDED  APPLICATION  1  RANGE  0-20,000 
LBS /ACRE 

?  0 

Crop  Residue  Applications 

You  are  asked  to  enter  the  number  of  organic  applications.   If 
there  were  no  applications,  enter  0  and  the  remainder  of  the 
input  is  bypassed.   If  organic  materials  were  applied,  the 
application  date,  depth  of  application,  and  content  are 
requested  for  each  application. 

THE  FOLLOWING  GUIDES  IN  ORGANIC  RESIDUE  DATA  INPUT 
INPUT  NUMBER  OF  ORGANIC  APPLICATIONS 

?  1 

ENTER  ORGANIC  RESIDUE  APPLICATION  DAY  (JULIAN)  FOR  APPLICATION 
NUMBER  1 

?  121 

ENTER  INCORPORATION  DEPTH  (CM)  OF  UNIFORM  ORGANIC  RESIDUE 
APPLICATION  1  (MUST  BE  A  POSITIVE  NUMBER)  ENTER  0.0  FOR  SURFACE 
APPLICATION 

?  6 
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ENTER  CARBON: NITROGEN  RATIO  OF  ORGANIC  MATTER  ADDED  IN 
APPLICATION  1  TYPICAL  RANGE  15-80  (C:N  RATIO  OF  40:1  IS  ENTERED 
AS  40.0) 

?  40.0 


ENTER  ORGANIC  NITROGEN  ADDED  (LBS/ACRE)  AS  OVEN  DRY  WEIGHT 
TYPICAL  RANGE  0-100  LBS/ACRE 

?  56 

The  three  previous  questions  are  repeated  for  each  organic 
application. 

ENTER  DEPTH  (CM)  OF  RESIDUE  ON  SURFACE 
TYPICAL  RANGE  0-10  CM 

?  5 

ENTER  FRESH  RESIDUE  REFLECTION  COEFFICIENT 

?  .8 

ENTER  RESIDUE  IDENTIFICATION  CODE 

0  IF  USER  WANTS  TO  ENTER  S ,  THE  FRACTION  OF  TOTAL 
RESIDUE  WHICH  IS  LABILE 

1  IF  CORN  STOVER 

2  IF  WHEAT  STRAW 

3  IF  CELLULOSE 

4  IF  SOYBEANS 

5  IF  S  IS  TO  BE  CALCULATED  FROM  THE  EQUATION 
S  =  0.07  +  l.ll(l./ACNl)**l./3. 

Note:   ACN1  is  the  carbon: nitrogen  ratio  of  residue  added 
(Data  Group  20) . 

When  the  organic  residue  data  input  has  been  completed,  the  data 
are  written  to  disk  file  and  program  control  is  returned  to  the 
operating  system. 
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RUNNING  THE  SOIL  TEMPERATURE  SUBMODEL 


This  section  assumes  that  the  user  has  successfully  installed  the 
soil  temperature  submodel  on  the  computer  and  has  been  able  to 
generate  the  required  input  data  files  AIR  and  INPUT1 .   Routine 
operation  of  the  submodel  involves  computer  runs  with  various 
climate  input  data  files  under  different  management  schemes  and 
initial  conditions.   The  user  is  responsible  for  the  file 
maintenance  and  storage  associated  with  the  many  computer  runs 
involved  in  a  complex  study.   The  file  naming  conventions  used 
in  this  model  require  the  user  to  rename  output  files  after  each 
run  and  to  make  certain  that  the  correct  input  files  are  named 
AIR  and  INPUT1 . 

Once  the  common  input  file,  INPUT1 ,  has  been  created,  related 
versions  with  minor  modifications  in  management  and  initial 
conditions  can  be  created  using  the  EDLIN  or  similar  text 
editor.   Details  concerning  specific  data  items  in  this  file  can 
be  found  in  the  "Common  Input  File  for  Soil  Temperature"  section 
of  this  manual.   Creation  of  additional  max-min  air  temperature 
files  can  be  done  the  most  efficiently  by  using  the  interactive 
data  setup  program. 

The  starting  day  for  the  soil  temperature  runs  should  correspond 
to  the  starting  day  used  in  the  general  NTRM  model.   For 
example,  if  the  soil  temperature  simulation  starts  on  day  121, 
then  the  general  NTRM  run  should  also  start  on  day  121. 

Both  the  soil  temperature  and  air  temperature  output  files  are 
structured  for  direct  input  into  the  general  NTRM  model.   Be 
sure  to  enter  the  correct  file  names  into  NTRFILE  before  running 
general  NTRM.    The  common  output  file,  OUTPUT,  is  saved  for 
listing  or  viewing  at  a  later  date .   Be  sure  to  rename  this  file 
if  you  want  to  save  the  results  from  each  run. 
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COMMON  INPUT  FILE  FOR  SOIL  TEMPERATURE  SUBMODEL 


Data  Group  1- - 
Starting  Day  and 
Number  of  Nodes 


The  soil  temperature  submodel  (Gupta  et  al .  1982)  requires  two 
input  files  and  produces  daily  mean  soil  temperature  and  hourly 
air  temperature  output  files  suitable  for  direct  input  into  the 
main  NTRM  model  (fig.  7).   Daily  maximum  and  minimum  air 
temperatures  must  be  supplied  in  a  file  structured  as  shown  in 
table  1.   In  addition,  a  general  input  file  (control  data), 
illustrated  in  figure  8,  must  be  prepared  as  follows: 


Column   Var . 

Var. 

Variable 

No .       type 

name 

description 

Line  1  of  2 

1-3      I 

JDS TART 

Julian  star 

run. 


Line  2  of  2 


1-3 


N 


Number  of  nodes. 


Data  Group  2 — 
Node  Depths  and 
Initial  Soil 
Temperatures 


Column   Var .    Var . 
No .       type    name 

Line  1  of  variable 

1-4  R  DP(1) 

5-10  R  TI 


11-14  R 

15-20  R 


DP(2) 

TJ 


Variable 
description 


Depth   of   first  node    (cm) . 

Initial    (at   time   zero  hour) 
soil    temperature    (°C) 
corresponding   DP(1) . 

Depth   of   second  node    (cm) . 

Initial    (at   time   zero  hour) 
soil    temperature    (°C) 
corresponding   to   DP(2) . 


Repeat   node    depth   and   soil    temperature    to 


^Model    is   currently   set    up   to   handle    38   nodes 
(see   example    in   section   "Running    the    Soil 
Temperature    Submodel"). 
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Column   Var .    Var. 
No.      type    name 


51-54 


55-60 


DP(6) 


TI 


Line  2  of  variable 


1-4 


5-10 


R 


DP(7) 


TI 


Variable 
description 

Depth  of  soil  temperature 
corresponding  to  DP(6) . 

Initial  (at  time  zero  hour) 
soil  temperature  (°C) 
corresponding  to  DP(6) . 


Depth  of  soil  temperature 
corresponding  to  DP(7). 

Initial  (at  time  zero  hour) 
soil  temperature  (CC) 
corresponding  to  DP(7). 


Repeat  depth  and  soil  temperature  to: 

51-54    R      DP(16) 

55-60    R      TI 

Line  3  of  variable 

Same  format  as  above  line- -continuation  of  same  data  up  to  DP(N) 
and  corresponding  TI . 
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SOIL 

TEMPERATURE 

SUBMODEL 


TO  GENERAL 
NTRM  MODEL 


Figure  7 

Input  and  output  configuration  for  Soil 

Temperature  submodel. 


GROUP  #6    THERMAL  DIFFUSIVITIES 


GROUP  #5   RESIDUE  AND/OR  CROP  COVER 
TREATMENTS 

GROUP  #4     CONSTANTS 


GROUP  #3     DAYS  IN  SIMULATION  AND 
HOURLY  PRINT  DAYS 


GROUP  #2     NODE  DEPTHS  AND  INITIAL 
SOIL  TEMPERATURES 


GROUP  #  1      STARTING  DAY  AND  NUMBER 
OF  NODES 


Figure  8 

Common  input  file  (control  data)  for  Soil 

Temperature  submodel. 
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Data  Group  3- -Days 
in  Simulation  and 
Hourly  Print  Days 


Table  1. 

File  structure  for  max-min  air  temperature  data  input  file. 


Record  number Type 


Contents 


1  Coded         Julian  date,  max  air 

temperature  (°F),  min  air 
temperature  (°F),  (spaced 
as  shown  below) . 

N  Coded        Repeat  for  each  day. 

Details  for  records  1  through  N: 

Begin  record  End  record 

AAA121     AAA62      AAA42 

1st  field      2nd  field      3rd  field 
6  characters    5  characters   5  characters 

Note:   Numbers  must  be  right  justified  in  these  fields. 


Column   Var . 


Var. 


No .       type    name 
Line  1  of  variable 


1-3 


KLMX 


4-6 


LPRINT 


Variable 
description 


Number  of  days  in  simulation. 
This  is  one  day  less  than  the 
number  of  days  in  the  max-min 
air  temperature  file. 
Example:  (Julian  stopping  day 
-Julian  starting  day 
(JDSTART))    =  KLMX.2 

If   LPRINT  =0,    the   program 

does   not   print  hourly    soil 

temperatures . 
If  LPRINT  =   1,    the   program 

prints   hourly   soil 

temperatures    for   every   day. 
If  LPRINT  =   2,    the   program 

prints  hourly   soil 

temperatures    for   a   limited 

number   of   days . 


^The  model   can  currently  handle   up   to    200 
days.      The   program  dimension   size   can  be   ex- 
panded  to   accomodate  more   days,    if   necessary. 
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Column   Var. 


No. 


type 


Var. 
name 


Variable 
description 


If  LPRINT  =  2,  the  following  lines  are  needed: 
Line  2  of  variable 


1-3 


NU 


Line  3  of  variable 
11-3      I       NDAY(l) 


Number  of  days  of  hourly 
print  (max  =  200) . 


Julian  dates  for  hourly 
print . 


Repeat  to: 

67-69     I       NDAY(23) 
Line  4  of  variable 


Same  format  as  above  card- -continuation  of  same  data 


Data  Group  4-  - 
Constants 


Data  Group  5-  - 

Residue  and/or  Crop   No 

Cover  Treatments 


Column 

Var. 

Var. 
name 

Variable 
description 

No. 

type 

Line  1 

of  1 

1-3 

R 

NX 

This  should 

be  1. 

Column 

Var. 

Var. 

Variable 

No. 

type 

name 

description 

Line  1 

of  1 

1-3 


R 


TR  EAT 


4-6 


MTH 


Defines  the  residue  and/or 

crop  cover  treatment. 
1  =  bare  soil 
3  =  50%  oat  residue  cover 
5  =  75%  oat  residue  cover 
7  =  99%  oat  residue  cover 
9  =  corn  +  50%  oat  cover 
11  =  corn  +  no  residue  cover 

This  should  be  1. 
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Data  Group  6-  - 

Column 

Var. 

Var. 

Thermal 

No. 

type 

name 

Diffusivities 

Line  1 

of  vari 

able 

1-6 

R 

ALP(l) 

7-12 

R 

(2) 

13-18 

R 

(3) 

19-24 

R 

ALP(4) 

25-30 

R 

(5) 

31-36 

R 

(6) 

Line  2 

of  vari 

able 

1-6       R      ALP ( 7 ) 


Repeat  to : 

73-78    R      ALP(19) 

Line  3  of  variable 


Variable 
description 


Thermal  diffusivity  (cmz/h) 


Thermal  diffusivity  (cm^/h) 


Same  format  as  above- -continuation  of  same  data  up  to  ALP(N) ,  N 
being  the  number  of  nodes  specified  in  Data  Group  1. 
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SOIL  TEMPERATURE  AND  AIR  TEMPERATURE 
FILES  (INPUT  TO  CENTRAL  NTRM  MODEL) 


These  input  files  are  normally  created  by  the  soil  temperature 
submodel.   They  are  described  here  in  detail  to  provide 
information  for  analysis  and  to  permit  inputs  from  other 
sources . 

Soil  temperature  file  (tape  20) 

Record    Type  Format  Variable  Contents 

1       Coded  (1X.F6.2,   DIST(l)  Depth  to  node  1  (cm). 

1       Coded  F6.2,       DIST(2)  Depth  to  node  2  (cm). 


1  Coded  F6.2,  DIST(20)  Depth  to  node  20  (cm) 
1  Coded  11X.F6.2,  DIST(21)  Depth  to  node  21  (cm) 
1       Coded     F6.2,       DIST(22)      Depth  to  node  22  (cm) 


Repeat 

to : 

1 

Coded 

2 

Coded 

3 

Coded 

F6.2)       DIST(30)      Depth  to  node  30  (cm). 

(13X.I3)    IDAY         Julian  day. 

(11X.F6.1,  TEMP(l)      Soil  temperature  for 

node  1  (°C) . 

Coded     F6.1,       TEMP(2)       Soil  temperature  for 

node  2  (°C) . 
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Repeat  to: 

3      Coded    F6 . 1 , 


TEMP(20)      Soil  temperature  for 
node  20  (°C) . 


Coded    11X.F6.1,   TEMP(21) 


Coded    F6.1, 


Soil  temperature  for 
node  21  (°C) . 


TEMP(22)      Soil  temperature  for 
node  22  (°C) . 


Repeat  to : 

3      Coded    F6.1) 


TEMP(30)      Soil  temperature  for 
node  30  (°C) . 


Record  1  is  read  once  per  run.   Records  2  and  3  are  repeated  in 
sequence  for  each  Julian  day. 


Air  temperature  file  (tape  21) 
Record   Type     Format     Variable 


1 

Coded 

(6X.I3) 

I  DAY 

2 

Coded 

(5X,F5.1, 

AIRT(l) 

2 

Coded 

F5.1, 

AIRT ( 2 ) 

2 

Coded 

F5.1, 

AIRT(3) 

Contents 

Julian  day. 

Hourly  air  temperature 
(0000-0100) . 

Hourly  air  temperature 
(0100-0200) . 

Hourly  air  temperature 
(0200-0300)  . 


Repeat  to: 

2      Coded    F5.1) 


AIRT (24)      Hourly  air  temperature 
(2300-0000)  . 


Sequence  of  two  records  repeated  for  each  Julian  day. 
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RUNNING  THE  CENTRAL  NTRM  MODEL 


This  section  assumes  that  the  user  has  sucessfully  installed  the 
general  NTRM  model  and  has  been  able  to  run  RNTRM2.EXE  and 
generate  the  input  files  required  to  run  this  model. 
The  general  configuration  of  NTRM  input  and  output  files  is 
shown  in  figure  9 .   Specific  details  of  the  file  structures  and 
data  formats  are  given  in  appendix  A. 

Names  for  the  input  and  output  files  must  be  specified  in  a  file 
called  NTRFILE.   The  model  uses  this  file  to  set  up  the  1-0  for 
a  particular  simulation  run.   The  user  should  rename  any  output 
files  which  must  be  saved  before  running  the  model  again. 
Otherwise,  the  file  will  be  overwritten  with  data  from  the  new 
run. 

The  soil  and  air  temperature  files  generated  by  the  soil 
temperature  submodel  can  be  used  as  direct  input  to  the  model. 
The  proper  file  names  need  only  be  specified  in  NTRFILE. 

The  control  switches  and  initial  soil  properties  are  contained 
in  the  common  input  file.   This  file  is  the  key  to  making 
various  types  of  runs  with  NTRM.   During  routine  applications, 
the  common  input  file  can  be  modified  with  the  text  editor  to 
establish  a  variety  of  possible  run  configurations  and  initial 
soil  properties .   The  user  should  become  familar  with  the 
locations  of  various  switches  and  control  parameters  in  this 
file. 

Output  from  the  model  occurs  in  three  basic  forms,  (1)  files 
intended  for  input  to  a  graphics  package,  (2)  tabular  general 
output  files  with  information  intended  for  a  line  printer,  and 
(3)  files  intended  for  input  to  other  simulation  models.   In 
general,  the  graphics  (plot)  files  can  be  examined  using  a  text 
editor  or  line  printer,  but  the  lack  of  detailed  formatting 
makes  this  process  more  difficult.   The  best  way  to  use  these 
files  is  via  a  plot  package  where  visual  comparisons  can  be  done 
quickly  and  easily. 

The  general  output  files  contain  information  such  as  soil 
constituent  masses  as  a  function  of  depth  and  time,  a  detailed 
mass  balance  summary  for  nitrogen,  details  of  each  water 
application  event,  top  and  root  stress  values  over  time,  and 
soil  water  contents  and  suctions  as  a  function  of  depth  and 
time . 
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FROM  SOIL  TEMP  SUBMODEL 


(SEE  APPENDIX  B) 
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AQUIFER 
SUBMODELS 


Figure    9 

Input  and  output  configuration  for  central 

NTRM  model. 


Whenever  a  new  soil  or  crop  configuration  is  established  in  the 
model,  the  model  results  should  be  checked  against  field 
validation  data  and  basic  experience.   Some  field  conditions  are 
difficult  to  simulate  with  only  a  casual  inspection  of  the  soil 
and  climate  properties.   You  could  easily  set  up  the  model  to 
simulate  a  set  of  conditions  considerably  different  from  those 
actually  occurring  in  the  field.   This  situation  occurs  because 
the  user  must  always  simplify  the  field  conditions  in  order  to 
fit  them  into  the  model.   It  is  easy  to  omit  certain  key  parts 
of  the  field  system,  while  others  are  defined  in  too  much 
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detail.   Experience  combined  with  a  reasonable  amount  of  field 
checking  are  required  in  order  to  set  up  a  valid  simulation. 

The  soil  horizon,  segment,  and  nodal  structure  used  in  the  NTRM 
model  are  illustrated  in  figure  10.   This  one -dimensional 
configuration  is  set  by  the  user  before  a  run  is  started.   The 
model  may  alter  the  spacing  in  the  plow  layer  depending  on  the 
initial  conditions  and  subsequent  inputs  to  the  system.   The 
basic  conventions  for  spacial  configurations  used  in  the  model 
are  outlined  below. 

Conventions : 

1.  Horizons  are  the  largest  subdivisions  (max=10) . 

2.  There  are  1  to  24  segments  per  horizon. 

--most  model  calculations  are  done  on  a  segment  basis, 
--there  can  be  a  maximum  of  24  segments  in  the  soil  column. 

3 .  Nodes  for  flow  calculations  occur  on  segment  boundaries  or 
are  evenly  spaced  between  boundaries  in  each  horizon. 
Nodal  spacing  is  constant  within  each  horizon. 

4.  The  user  specifies  the  number  of  horizons  and  the  depth  to 
the  lower  boundary  of  each. 

5 .  The  user  specifies  the  number  of  uniformly  spaced  segments 
per  horizon. 

6.  The  spacing  between  segments  is  known  as  DELX.   The  spacing 
between  nodes  is  known  as  DELXF  (in  most  cases,  DELX  = 
DELXF) . 

7 .  The  user  specifies  the  number  of  nodes  evenly  spaced 
between  horizon  boundaries . 

8.  Segment  (1)  represents  the  soil  surface  and  contains  any 
crop  residues  or  surface  applications  of  fertilizer. 


User  Guide  to  Horizons ,  Segments ,  and  Nodes 

1.  First  set  up  the  horizon  locations  based  on  available  soils 
data . 

2 .  Then  subdivide  the  horizons  into  segments  based  on  the 
detail  desired  therein. 

3.  Nodes  must  occur  on  each  segment  and  horizon  boundary. 

4.  Additional  nodes  can  be  added  between  horizon  boundaries 
if  necessary. 
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Soil  Surface 
Horizon  (1) 


Horizon  (2) 


DELX 


,d> 


Segment  (1) 
Segment  (2) 


Segment  (3) 


Segment  (4) 


Segment  (5) 


Segment  (6) 


Segment  (7) 


Horizon  (3) 

Segment  (8) 

, « i 

•        DELX  (9) 

Segment  (9) 

Horizon  (4) 

A 

• 

• 

•       DELX  (10) 

Segment  (10) 

S 
i ._ 

Nodes 

Figure    10 

NTRM  model   horizon, 

nodal    structure. 

seg 

ment ,    and 
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NTRM  MODEL  PLOTTED  OUTPUT 


A  considerable  volume  of  model  output  is  condensed  in  the  form 
of  two-  and  three-dimensional  plot  output  files.   The  model  is 
capable  of  producing  a  variety  of  plot  types  depending  on  user 
specifications  and  requirements.   A  summary  of  the  plots  that 
can  be  generated  from  each  plot  file  is  given  in  table  2 . 


Table  2. 

NTRM  model  plot  output  summary 


Plot  file: 


Plots  generated 


Tape  2  (Deep  Percolation  Plot  File) 


Tape  4  (Top  Growth  Plot  File) 


1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 
10. 
11. 
12. 


Cumulative  leachate  volumes  (cm) 

vs.  time. 

Incremental  leachate  volumes 

(cm)  vs .  time . 

NO3-N  cone,  in  leachate 

(mg/L)  vs.  time. 

NHJ-N  cone,  in  leachate 

(mg/L)  vs.  time. 

UREA-N  cone,  in  leachate 

(mg/L)  vs .  time . 

in  leachate 

time  . 
in  leachate 

time  . 

Mg"1-*"  cone .  in  leachate 
(mg/L)  vs .  time . 
HCO3  cone,  in  leachate 
(mg/L)  vs.  time. 
CI"  cone,  in  leachate 
(mg/L)  vs.  time. 
CO3  cone .  in  leachate 
(mg/L)  vs.  time. 

so: 


Ca   cone . 
(mg/L)  vs . 
Na   cone . 
(mg/L)  vs. 


4  cone .  in  leachate 
(mg/L)  vs.  time. 


1.  Total  dry  matter  (tops)  vs. 
time . 

2.  Leaf  area  (dm^/plant)  vs.  time. 

3.  Corn  grain  production  (kg/ha) 
vs .  time . 

4.  N  uptake  demand  (kg/ha)  vs. 
time . 


Tape  6  (Root  Growth  Plot  File) 


Total  net  root  length  produced 
(cm/cm^  soil)  vs.  time. 
Total  root  length  produced 
(era/cm^  soil)  vs.  time. 
Total  root  length  sloughed 
(cm/cm^  soil)  vs.  time. 
Total  root  mass  produced  (kg/ha) 
vs.  time. 
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Table  2- -Continued 

NTRM  model  plot  output  summary 


Plot  file: 


Plots  generated 


Tape  7  (Root  Distribution  Plot  File) 


Tape  12  (Soil  Nitrogen  Plot  File) 


Tape  13  (3-D  Plot  File) 


1.   Crop  root  distribution  vs. 
depth  and  time. 

1.  Total  NHJ-N+NO3-N  in  soil 
(kg/ha)  vs.  time. 

2.  Daily  plant  uptake  of  N 
(kg/ha/day)  vs.  time. 

3.  Daily  N  leached  (kg/ha/day) 
vs.  time. 

4.  Daily  N  denitrified 
(kg/ha/day)  vs.  time. 

5.  Daily  N  fixed  (kg/ha)  vs. 
time . 

6.  NO3-N  on  soil  surface 
(kg/ha)  vs.  time. 

7.  NHJ-N  on  soil  surface 
(kg/ha)  vs.  time. 

1.  Root  length  (cm/cm  soil) 
vs.  depth  (cm)  and  time. 

2.  Soil  temperature  (°C)  vs. 
depth  (cm)  and  time. 

3.  Soil  water  content 

3     T 

(cm  /cm  )  vs.  depth  (cm) 
and  time. 

4.  Soil  suction  (cm  H2O)  vs. 
depth  (cm)  and  time. 

5.  Soil  bulk  density  (g/cm-5) 
vs.  depth  (cm)  and  time. 

6.  NO3-N  cone.  (mg/L) 

vs.  depth  (cm)  and  time. 

7.  NH^-N  cone.  (mg/L)  vs. 
depth  (cm)  and  time. 

8.  UREA-N  cone.  (mg/L) 

vs.  depth  (cm)  and  time. 

9.  Ca++  cone.  (mg/L)  vs. 
depth  (cm)  and  time. 

10.  Na+  cone.  (mg/L)  vs.  depth 
(cm)  and  time. 

11.  Mg++  cone.  (mg/L)  vs.  depth 
(cm)  and  time. 

12.  HCO3  cone.  (mg/L)  vs.  depth 
(cm)  and  time. 

13.  CI"  cone.  (mg/L)  vs.  depth 
(cm)  and  time. 

14.  CO3  cone.  (mg/L)  vs.  depth 
(cm)  and  time. 

15.  SO4  cone.  (mg/L)  vs.  depth 
(cm)  and  time. 


' 
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Table  2- -Continued 

NTRM  model  plot  output  summary 


Plot  file: 


Plots  generated 


Tape  14  (Soil  Surface  Plot  File) 


Tape  15  (N15  Plot  File) 


1. 
2. 

3. 

4. 
5. 

6. 

7. 

1. 
2 . 
3. 
4. 
5. 
6. 

7. 


Random  roughness  (cm)  vs.  time. 
Bulk  density  of  uppermost 
segment  (g/cm  )  vs.  time. 
Surface  residue  (mass) 
(t/ha)  vs.  time. 
Percent  cover  vs.  time. 
Total  water  content  (cm)  of 
canopy  vs .  time . 
Total  water  content  (cm)  of 
residue  (mulch)  vs.  time. 
Surface  depression  storage 
(cm)  vs.  time. 

NO3-N  (%  N15)  vs.  time. 
SNHt-N  (%  N15)  vs.  time. 
RNHJ-N  (%  N15)  vs.  time. 
Soil  OM  (%  N15)  vs.  time. 
Residues  (%  N15)  vs.  time. 
Plant  N  uptake  (%  N15)  vs. 
time . 

time . 

vs . 


time . 

Leached  N  (%  N15)  vs. 
Denitrified  N  (%  N15) 
time . 


Tape  17  (Evapotranspiration  Plot  File) 


Daily  potential  evaporation 
from  soil  surface  (cm/mo) 
vs.  time. 
Daily  potential 
transpiration  (cm/mo)  vs. 
time . 

Daily  potential  ET 
(cm/mo)  vs.  time. 
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COMMON  INPUT  FILE  FOR  CENTRAL  NTRM  MODEL 


The  basic  data  groups  in  the  common  input  file  needed  to  run  the 
general  NTRM  model  are  illustrated  in  figure  11.   Certain  groups 
(see  separate  group  descriptions)  must  be  omitted  when  the 
nitrogen  portion  of  the  model  is  bypassed  (NBYPAS  =  1) .   Group 
11  must  be  included  when  the  salt  portion  of  the  model  is 
bypassed  (SBYPAS  =1).   At  all  others  times,  group  11  is 
omitted. 


In  a  similar  manner,  group  12  is  included  only  when  flow  data 
are  to  be  read  from  cards  (ITEST  =  1).   When  ITEST  =  0,  flow 
data  are  read  from  tape  1 . 

Group  14,  the  restart  data,  is  not  currently  used  in  NTRM.   The 
feature  may  be  reactivated  at  a  future  date. 


Group  21,  Root  Distribution,  is  present  only  when  IVARR  =0.   If 
the  root  growth  model  is  used,  IVARR  =  2,  group  21  is  omitted. 

Group  22,  Plant  Uptake  of  Nitrogen,  is  present  only  when  CR0P=3 . 


Group  23,  Crop  Growth  Information  (corn)  is  present  only  when 
CROP  =  4. 


Data  groups  24  through  26  are  the  input  data  for  the  unsaturated 
flow  submodel. 

These  data  groups  are  similar  to  those  used  by  previous 
independent  versions  of  Unsaturated  Flow.   Modifications  have 
been  made  in  certain  groups  to  accommodate  the  layered  soil  and 
modified  boundary  conditions . 


Data  Group  1- -Title   Column   Var . 

No .       type 


Var 


name 


Variable 
description 


Line  1  of  1 


■  ■ 


1-80 


TITLE 


Any  title  desired  on  printed  output 
This  character  string  will  also 
appear  in  right  margin  of  plotted 
output. 


Data  Group  2- - 
Control  Data 


Column   Var . 


Var 


No. 

type 

name 

Line 

1 

of  4 

1-5 

X 

none 

6-10 

X 

none 

Variable 
description 

Not  used  by  program. 
Not  used  by  program, 
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■ 


GROUP  #26    SOIL  PHYSICAL  PROPERTIES 


K~ 


GROUP  #25     10  OPTIONS  AND  RUN 

CONDITIONS  FOR  UNSATURATED 
FLOW 


GROUP  #24    SOIL  SURFACE  AND  LOCATION 
DATA 


GROUP  #23    CROP  GROWTH  INFORMATION 


GROUP #22    PLANT  UPTAKE  OF  NITROGEN 

GROUP  #21     ROOT  DISTRIBUTION 

GROUP  #20    ORGANIC  RESIDUE 
APPLICATIONS 


GROUP  #19     FERTILIZER  APPLICATIONS 


GROUP  #18    TILLAGE  UPDATES 


GROUP  #17    TILLAGE  DATES,  DEPTHS 
AND  CODES 


GROUP  #16    ORGANIC  RESIDUE 

APPLICATION  DATES 


GROUP  #15     FERTILIZER  APPLICATION 
DATES 


c 


GROUP  #14     RESTART  DATA 


GROUP  #13     INITIAL  SOIL  ANALYSIS 


GROUP #12     FLOW  DATA 


GROUP  #11     SALT  BY-PASS  DATA 


GROUP  #10     RAIN  WATER  ANALYSIS 


GROUP  #9     IRRIGATION  WATER  ANALYSIS 


GROUP #8    TEMPERATURE  DATA 


GROUP  #7    COEFFICIENTS  FOR  TRANSITION 

STATE  NITRIFICATION  EQUATIONS 


GROUP #6    NODAL SPACINGS PER  SEGMENT 


GROUP  #5    SEGMENTS  PER  HORIZON 


GROUP  #4    SOIL  HORIZON  DEPTHS 


L 


GROUP  #3  TEMPERATURE  HORIZON  DEPTHS 


GROUP #2  CONTROL  CARDS 


GROUP  #1  TITLE  CARD 


OPTIONAL 


REQUIRED 

Figure  11 

Basic  data  groups  in  the  common  input 

file  needed  to  run  the  central  NTRM  model. 
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Column 

Var. 

Var. 

Variable 

No. 

type 

name 

description 

11-15 

R 

RSN 

Rate  of  bre< 

16-20 


21-25 


R 


R 


CH 


CHI 


26-30 

X 

none 

31-35 

X 

none 

36-40 

X 

none 

41-45 

X 

none 

46-50 

X 

none 

51-55 

R 

FACT 

56-60 


R 


SDATE 


(lb/acre/yr) .   (Shaffer  nitrogen 
transformation  (NT)  model  only.) 

Convergence  criterion  for  nitrogen 
transformation  routine.   See  note 
for  details .   Suggested  values  are 
in  the  range  0.1-10.0  ppm  nitrogen. 
(Shaffer  NT  model  only) 

Convergence  criterion  for  chemical 

9 

equilibria  subroutine.    Suggested 
values  are  in  the  range  1*10 
-1-10"6  moles/L. 

Not  used  by  program. 

Not  used  by  program. 

Not  used  by  program. 

Not  used  by  program. 

Not  used  by  program. 

Plant  N-uptake  constant  (/zg/cm)  for 
use  with  consumptive -use  values  (if 
used) . 

Emergence  date  (1-365). 


Line  2  of  4 


1-5      I 
6-10     I 


LL 

MM 


Julian  starting  day  of  run  (1-365) 

Julian  termination  day  of  run 
(1-365) . 


^The  constant  CH  determines  the  accuracy 
of  the  nitrogen  changes  computed  in  sub- 
routine TRNSFM  (Shaffer  model).   A  Smaller 
value  will  increase  accuracy  and  program 
execution  time. 

^The  constant  CHI  determines  the  accuracy 
of  mass  changes  in  the  system  computed  in 
subroutine  XCHANGE.   A  Smaller  value  will 
increase  accuracy  and  program  execution 
time. 
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Column 

Var. 

Var. 

Variable 

No. 

type 

name 

description 

11-15 

I 

0 

Number  of  s< 

be  in  range  1  to  10  inclusive)  (see 
figure  5) . 

16-20     I       CROP         If  =  1,  transpiration  used  to 

determine  N-uptake  by  plants. 

If  =  2 ,  reserved  for  future  use. 

If  =  3,  read  plant  N-uptake 

data  (see  data  groups  16  and  17). 

If  =  4,  crop  growth  submodel 
used  (for  corn) . 

21-25      I       TO         Number  of  temperature  horizons 

(must  be  in  range  1  to  4  inclusive) . 
(Leave  blank  if  ITEMP  =  1,  see  group 
2.) 

26-30      I       NT         Number  of  weekly  temperature  data 

lines  (must  =  0  for  NBYPAS  =  1,  see 
this  group  and  group  8) .  (Leave 
blank  if  ITEMP  =  11,  see  group  2.) 

31-35      I       IPRINT      Print  interval  (days)  for  predicted 

soil  chemistry  data. 

36-40     I      JPRINT     Time  interval  (1-10)  within  which 

the  IPRINT  and  IMASS  print  options 
are  to  be  activated. 

41-45     I      INK        Write  interval  (days)  for  output  to 

tape  2  (see  "Pan  Evaporation,  Wind 
and  Solar  Radiation  File"). 

Not  used  by  program. 

Not  used  by  program. 

Not  used  by  program. 

If  =  0,  waterflow  data  are 
read  from  tape  1  or  computed 
internally. 

If  =  1,  steady  state  waterflow 
data  are  read  (see  group  12). 


46-50 

X 

none 

51-55 

X 

none 

56-60 

X 

none 

61-65 

I 

ITES' 

U'6 


Column 

Var. 

Var. 

Variable 

No. 

type 

name 

description 

66-70 

I 

I  CAP 

If  =  1,  upw< 

allowed  in  plug  flow  submodel. 
If  =  0,  upward  movement  of  water  not 
allowed  in  plug  flow  submodel. 

71-75      I       NCAP       Number  of  nodes  associated  with 

capillary  fringe  (from  the  bottom 
node  upward) . 

76-80      I       KYEAR      Year  for  this  run  (for  example,  if 

this  is  the  third  year  being  run, 
KYEAR  =  3 ) . 


Line  3  of  4 

1-5       I       KSTOP 


Year  when  run  terminates  (for 
example,  if  run  Is  to  terminate 
after  33  years  and  starts  with  year 
3,  KSTOP  =  5) . 


6-10      I       IMASS        Print  interval  (days)  for  summary  of 

nitrogen  balance  for  system. 

11-15     I      IPRIN8      If  =  0,  printing  of  below  is 

suppressed. 
If  =  1,  basic  control  data  are 
printed. 

16-20     I      IPRIN7      If  =  0,  printing  of  below  is 

suppressed. 
If  =  1,  input  data  are  printed. 

21-25     I      NBYPAS      If  =  0,  nitrogen  transformations  are 

considered  (Shaffer  model  or 
Molina  model) . 
If  —  1,  nitrogen  transformations  are 
bypassed,  salt  reactions  remain  in 
model  if  SBYPAS  =  0. 


26-30    I 


SBYPAS 


31-35     I 


36-40    I 


NROTA 


NCLI 


If  =  0,  salt  reactions  are 
considered. 

If  =  1,  salt  reactions  are  bypassed; 
nitrogen  transformations  remain 
in  the  model  if  NBYPAS  =  0. 

Number  of  years  in  crop  rotation 
cycle . 

Number  of  years  in  climate  cycle. 


4  9 


Column 

Var. 

Var. 

Variable 

No. 

tvpe 

name 

description 

41-45 

I 

INFIL2 

Starting  fi' 

46-50     I 


51-55     I 


56-60     I 


INFIL1 


IC0NT1 


IWRI3 


(used  when  making  a  restart  on 
tape  2  containing  previous  files) . 

Starting  file  on  input  tape  1 
(used  when  more  than  1  year's  data 
are  to  be  input  to  flow 
calculations) . 

If  =  1,  a  new  set  of  water 
application  dates  is  read  in 
each  year.   Also,  no  records 
can  be  skipped  on  tape  1. 

If  =  0,  entire  3-D  plot  file 

(tape  13)  is  written. 
If  =  1,  first  part  of  3-D  plot 

file  is  written. 


61-65     I 


[OP 


66-70    I 


IOPN 


71-75     I 


IPC02 


76-80     I 


IREK 


If  =  1,  lime  precipitation  is 

checked  in  relation  to  available 
pore  space . 

If  =  0,  option  is  bypassed. 

If  =  1,  transition-state 

nitrification  subroutine  used 
(Shaffer  and  Molina  NT  models) . 
Coefficients  must  be  read  (see 
group  7) . 

If  =  0,  statistical  nitrification 
equation  used  (Shaffer  model 
only) . 

If  =  1,  user  supplies  soil  pC02 

values  (see  group  13) . 
If  =  0,  model  estimates  soil  pC02 

values . 

If  =  1,  cation  exchange  coefficients 
for  Ca-Na  and  Ca-Mg  ion  exchange 
are  read  from  soils  analysis  data. 

If  =  0,  default  values  within 

program  are  used.   (Ca-Na  =  0.707 
(moles/L)"1/2  and  Ca-Mg  = 
0.670  (unitless  ratio)). 
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Column   Var . 
No .       type 

Line  4  of  4 


Var, 


name 


Variable 
description 


1-5 


IVARR 


6-10 


11-15 


16-20    I 


21-25     I 


ITEMP 


IDEMD 


I  FLOW 


INIT 


Determines  how  root  growth  will  be 

simulated. 
If  =  0,  constant  root  distribution 

from  data  group  21. 
If  =  2,  root  growth  model  of 

Shaffer  and  Clapp  (1982)  is  used. 

Soil  temperature  control. 

If  =  0,  weekly  soil  temperatures 
are  read  in  (group  8) . 

If  =  1,  soil  temperature  model  of 
Gupta  et  al.  (1982)  is  used.   (See 
section  on  soil  temperature 
submodel  for  instructions . ) 

Demand  control  for  crop  growth 

model . 
If  =  0,  crop  growth  is 

interactive  with  nitrogen 

availability. 
If  =  1,  crop  growth  is  based  on 

ideal  curve  with  respect  to 

nitrogen. 

Unsaturated  flow  model  control. 
If  =  0,  unsaturated  flow  data  must 

be  supplied  via  tape  1  or  data 

group  12 . 
If  =  1,  unsaturated  flow  program 

in  model  is  used. 

Nitrogen  transformation  model 

control . 
If  =  0,  Shaffer  NT  model  used. 
If  =  1,  Molina  NT  model  used. 


26-30     I 


JSET 


Root  growth  interaction  control 

switch  (Shaffer  root  growth  model 

only) . 
If  =  0,  root  growth  interactive  with 

soil  conditions . 
If  =  1,  root  growth  is  ideal  with 

respect  to  soil  conditions. 
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Column   Var . 
No .       type 

31-35     I 


36-40 


Var. 


name 


N15 


IEVAP 


Variable 
description 

N15  nitrogen  write-control  switch. 
If  =  0,  N15  data  not  output  to 

tape  15 . 
If  =  1,  N15  data  written  to  tape  15 

on  a  daily  basis. 

If  =  0,  potential  evaporation  set 

equal  to  zero . 
If  =  1,  potential  evaporation 

estimated  from  pan  evaporation 

and  ALAD. 


41-45     I      IDAYP8      If  =  0,  semimonthly  pan  evaporation 

read  from  tape  81. 
If  =  1,  daily  pan  evaporation,  wind, 
and  solar  radiation  read  from  tape 
81. 


46-50    I 


ITMOD 


If  =  0,  tillage  data  not  read. 
If  =  1,  tillage  data  read  from 
group  18 . 


51-55     I 


ITMSUB 


56-60    X      none 
61-65     I       IMETH 


66-70    I 


71-75     I 


1  PLUG 


JEVAP 


If  =  0,  tillage  data  read  from 
group  18  used  in  tillage 
calculations . 

If  =  1,  tillage  submodels  used  to 
make  tillage  calculations . 

Not  used  by  program. 

If  =  0,  independent  and  cross- 
product  terms  used  in  root 
submodel . 

If  =  1,  multiplicative  terms  used. 

If  =  2,  minimum  function  used. 

If  =  0,  unsaturated  flow  model  is 

used. 
If  =  1,  plug  flow  model  is  used. 

Number  of  nodes  for  surface 
evaporation  (plug  flow  model  only) 


76-80     I 


IWRI2 


If  =  0,  deep  percolation  data  are 

written  in  binary  form. 
If  =  1,  deep  percolation  data  are 

written  in  coded  form. 
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Data  Group  3- - 
Temperature  Horizon 
Depths 


(Omit  if  NBYPAS  =  1 ,  or  ITEMP  =  1;  see  Data  Group  2.) 


Data  Group  4- -Soil 
Horizon  Depths 


1 


Data  Group  5 — 
Segments  per 
Horizon 


Column   Var . 
No .       type 

Line  1  of  1 


Var. 
name 


1-5 


R 


6-103    R 


THOR(J) 
THOR(J) 


(Must  always  be  present.) 

Column   Var .    Var . 
No .       type 

Line  1  of  1 


1-5 


R 


4 


6-10 


Column   Var. 


No 


type 


name 

HOR(J) 
HOR(J) 

Var. 
name 


Line  1  of  1 


Var i  able 
description 


Depth  (cm)  from  surface  to  bottom 
of  first  temperature  horizon. 

Depth  (cm)  from  surface  to  bottom 
of  second  temperature  horizon  (if 
it  exists) .   (J  =  subscript  for 
each  horizon. ) 


Variable 
description 


Depth  (cm)  from  surface  to  bottom 
bottom  of  first  soil  horizon. 

Depth  (cm)  from  surface  to  bottom 
of  second  soil  horizon.    (J  = 
subscript  for  each  horizon.) 

Variable 
description 


1-5 
6-10 


ISEGN(J)     Number  of  segments  in  Jth 

horizon. 
ISEGN(J+1)   Number  of  segments  in  horizon 

J+l. 


^The  program  allows  the  inclusion  of  up  to 
four  temperature  horizons.   If  included,  the 
depths  for  numbers  3  and  4  would  be  in 
columns  11-15  and  16-20,  respectively. 
^The  program  allows  inclusion  of  up  to  10 
soil  horizons.   If  included,  the  depths  for 
numbers  3,  4,  5,  6,  7,  8,  9,  and  10  would  be 
in  columns  11-15,  16-20,  21-25,  26-30,  31-35, 
36-40,  41-45,  46-50,  respectively. 

^If  only  one  soil  horizon  is  desired,  only 
one  need  be  included. 


5  3 


Data  Group  6 --Nodal   Column   Var .    Var .         Variable 
Spacings  per  Segment  No.       type   name        description 

Line  1  of  1 

1-5       I       INODE(J)     Number  of  nodal  spacings  per 

segment  in  Jth  horizon. 

6-10      I       IN0DE(J+1)   Number  of  nodal  spacings  per 

segment  in  J+l  horizon. 

Data  Group  7--        (Needed  if  IOPN  =  1,  see  Data  Group  2.) 

Coefficients  for 

Transition  State      Column   Var .    Var.        Variable 

Nitrification         No.       type    name         description 

Equations 

Line  1  of  1 

1-5      R      AINTN       Y  intercept  of  activation  energy 

vs.  ionic  strength  relationship 
(default  =44.6  kcal/mole) . 

6-10      R       SLOPEN       Slope  of  activation  energy 

vs .  ionic  strength  relationship 
(default  =12.3  kcal-L/mole2) . 

Data  Group  8--        (Omit  if  NBYPAS  =  1,  or  ITEMP  =  1;  see  Data  Group  2.) 

Weekly  Soil 

Temperature  Date"     Column   Var .    Var.        Variable 

No.       type   name        description 

Line  1  of  NT7 

1-2      X      none        Card  identification  (not  read 

by  program) . 

3-10  R  TT(K)  Average  weekly   temperature    (°C) 

for   uppermost    (first)    temperature 
horizon. 


"The   first   temperature  data   should    correspond 

to   the   week  beginning  with  the   start   of   the 

run. 

'NT  is  the  number  of  temperature  data  lines. 

Each  line  contains  weekly  data  for  all 

temperature  horizons. 
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Column   Var. 


No. 
11-20* 


type 


Var. 
name 


TT(K) 


Variable 
description 

Average  weekly  temperature  (°C) 
for  second  temperature  horizon. 
(K  =  subscript  for  weekly 
temperatures . ) 


Data  Group  9-- 
Irrigation  Water 
Analysis 


Column 

Var. 

Var. 

Variable 

No. 

type 

name 

description 

Line  1 

of  1 

1-5 

R 

6-10 

R 

11-15 

R 

16-20 

R 

21-25 

R 

26-30 

R 

31-35 

R 

36-40 

R 

41-45 

R 

ANH3(1) 

AN03(1) 

CA(1) 

ANA(l) 

AMG(l) 

HC03(1) 

CL(1) 

C03(l) 

S04(l) 


(all  units  in  meq/liter  of 
the  species  shown) 


nhJ 


N03 


Ca" 


+ 


Na 

Mg++ 

HCO3 

CI" 

CO3 

SO4 


Note:   All  irrigation  water  is  assumed  to  have  the  same 
analysis . 


Data  Group  10- -Rain 
Water  Analysis 


Column   Var . 
No .       type 

Line  1  of  1 


1-5 


R 


Var. 
name 


ANH3(1) 


Variable 
description 


(all  units  in  meq/liter  of 
the  species  shown) 

NHj 


°The  program  allows  inclusion  of  up  to  four 
temperature  horizons.   If  included,  tempera- 
ture data  for  horizons  3  and  4  would  be  in 
columns  21-30  and  31-40,  respectively. 
'If  only  one  temperature  horizon  is  desired, 
only  one  need  be  included. 
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Column   Var.    Var.        Variable 


No. 

type 

name 

description 

6-10 

R 

AN03(1) 

NOj 

11-15 

R 

CA(1) 

Ca++ 

16-20 

R 

ANA(l) 

Na+ 

21-25 

R 

AMG(l) 

Mg++ 

26-30 

R 

HC03(1) 

HCO3 

31-35 

R 

CL(1) 

CI" 

36-40 

R 

C03(l) 

CO3 

41-45 

R 

S04(l) 

SO4 

Note:   All  rain  water  is  assumed  to  have  the  same 
analysis . 

Data  Group  11- -Salt   (Needed  if  SBYPAS  =  1  and  NBYPAS  =  0;  see  Data  Group  2.) 
Bypass  Data 

Column  Var .    Var.        Variable 

No.  type   name        description 

Line  1  of  Q-l 

1-5       R       SERAT8(J)    Fraction  of  total  NH4  that  is 

soluble   NHj. 

6-10  R  U(J)  Ionic    strength    (moles/L) 

--needed  only  when   I0PN=   1. 

11-15  R  ACTCA(J)  Ca++   activity    (moles/L) 

--needed  only  when   IOPN  =   1. 

Data   Group    12—  (Needed  only   if   ITEST  =   1;    if   ITEST  =   0,    tape    1  must  be   used 

Waterflow  Data  (IFLOW  =   0;    see   Data  Group    2)    or   internal   unsaturated   flow  model 

can  be    activated    (IFLOW  =   1;    see   Data  Group    2.) 


Column   Var. 

Var. 

Variable 

No.       type 

name 

description 

Line  1  of  NS10 

1-10  R  CMH201(J)         Volume    of  water    in   segment    (cm3) 


l^NS   is    the  number   of   waterflow  data   lines. 
There  must    be   a  line   for   each  segment 
including   No.    1. 
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Column   Var.    Var. 
No.      type 

11-20    R 


21-30    R 


name 


MOISIN(J) 


M0IS08(J) 


Variable 
description 

Water  movement  across  upper 
boundary  of  segment  (cm  /time 
step)  .  H 

Water  movement  across  lower  boundary 

o  11 

of  segment  (cm  /time  step). 


31-40     R      TEN(J)       Average  water  tension  (cm)  over  time 

step  (tensions  shown  as  negative) . 


41-50  R 


U(J) 


Plant  water   uptake    (cm/time    step 
/soil    segment) . 


Data   Group    13 
Initial    Soil 
Analysis 


Column        Var 


No 


Line   1   of  NC 


12 


1-5 


Var. 


type        name 


ANH3(L) 


Variable 
description 


(Units    in  meq/liter    of 
species    shown  unless 
otherwise   noted.) 

NHj    (soluble    form) . 
WARNING!!      Never   use    total 
NH.£  here. 


6-10 

R 

AN03(L) 

NO3 

11-15 

R 

UREA(L) 

UREA 

16-20 

R 

CA(1) 

Ca++ 

21-25 

R 

ANA(l) 

Na+ 

26-30 

R 

AMG(l) 

Mg++ 

31-35 

R 

HC03(1) 

HCO3 

36-40 

R 

CL(1) 

CI" 

41-45 

R 

C03(l) 

CO3 

46-50 

R 

S04(l) 

SO4 

^Positive    sign   for   downward    flow,    negative 

for  upward   flow. 

l^NC   is    the   number   of    soil   horizons. 
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Column 
No. 

51-55 


76-80 


Var. 
type 

R 


66-70    R 
71-75    R 


R 


Var. 
name 

EC(1) 


56-60    R      XX5(1) 
61-65     R      CAL(l) 


BD(1) 

SAMT(l) 

CN1(1) 


Variable 
description 

Exchange  capacity  (meq/100  g 
soil) . 

Gypsum  (meq/100  g  soil) . 

If  =  0.0,  no  lime  is  present  in 

soil . 
If  =  1.0,  lime  is  present  in  soil. 

o 

Bulk  density  (g/cm   of  soil) . 

Soil  organic  nitrogen  (/£g/g 
soil).   (Shaffer  NT  model  only.) 

Carbon-nitrogen  ratio  of  soil 
organic  matter  (for  example,  40:1 
would  be  entered  as  40.0). 
(Shaffer  NT  model  only.) 


Line  2  of  set 


1-5       R      XTRCT(l) 
6-15      R      PC02(1) 


16-20     R 


21-25     R 


26-30     R 


31-35     R 


AKCS(l) 


AKCM(l) 


CONCl(l) 


C0NC2(1) 


g  water/g  soil  in  soil  extract. 

CO2  partial  pressure  in  atmospheres 
(needed  only  if  IPC02  =  1). 

Ca-Na  exchange  coefficient  (needed 
only  if  IREK  =  1)  default  =  0.707 
(moles/L)"1/2. 

Ca-Mg  exchange  coefficient  (needed 
only  if  IREK  =  1)  default  =  0.670 
(unitless  ratio) . 

Microbial  mass  (/ig/g  soil,  as 
carbon).   (Molina  NT  model  only.) 

Nutrient  humus  (/ig/g  soil,  as 
carbon).   (Molina  NT  model  only.) 


Data  Group  14- 
Restart  Data 


(The  current  version  of  NTRM  does  not  allow  a  restart; 
this  option  may  be  reactivated  later.) 

This  data  group  is  written  by  the  program  during  a  previous  run. 
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Data  Group  15- - 

(Must  always  be 

present . ) 

Fertilizer 

Application  Dates 

Column   Var . 

Var. 

No .       type 

name 

Line  1  of  l13 

1-5       R 

ITOT 

6-10      I 

IADD(K) 

Data  Group  16- - 
Organic  Residue 
Application  Dates 


Data  Group  17-- 
Tillage  Dates, 
Depths,  and  Codes 


11-1514   I 


(Omit  if  NBYPAS 

Column   Var . 
No .       type 

Line  1  of  1 

1-5      I 

6-10     I 


11-1515   I 


Variable 
description 


Number  of  fertilizer  applications. 

IADD(K)      Julian  date  of  first  application. 

IADD(K)      Julian  date  of  second  application. 
(K  =  subscript  for  fertilizer 
applications . ) 


1 ;  see  Data  Group  2 . ) 


Var. 


name 


J  TOT 


IORNAP(K) 
IORNAP(K) 


Column        Var . 
No .  type 

Line    1    of    1 

1-5  I 


Var. 
name 


ITIL 


Variable 
description 


Number   of   organic   nitrogen 
(residue)    applications. 

Julian  date   of   first   organic 
nitrogen  application. 

Julian   date    of   second   organic 
nitrogen  application. 
(K  =   subscript   for   nitrogen 
applications . ) 

Variable 
description 


Number   of   tillage   events    (max  =6) 


l^A  second    line  may  be   used   if  needed    (that 
is,    if   the   number   of   applications   exceeds    15, 
in  which  case   the    15    format    still   applies). 

l^The   program  allows   inclusion   of   up   to    25 
fertilizer   application  dates.      The    same 
format,    15,    is   used    for   any   additional   dates 
desired. 

l^This   program  allows   inclusion  of   up   to    5 
organic   nitrogen  application  dates.      If 
included,    application  dates   for   applications 
3,    A,    and   5  would   be   in   columns    16-20,    21-25, 
and   26-30,    respectively. 
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Column   Var.    Var. 
No. 

6-10 

11-16 


17 


type   name 

I       ITILL(K) 
R       TILLD(K) 


Variable 
description 

Julian  date  of  Kth  tillage  event. 

Depth  (cm)  from  surface  to  bottom 
of  tillage. 


I       ITCODE(K) 
Repeat  sequence  up  to  6  times . 


Tillage  code  (1  =  moldboard  plow, 
2  =  chisel  plow,  3  =  disk). 


Data  Group  18-- 
Tillage  Updates 


(Omit  if  ITMOD  ==  0;  see  Data  Group  2.) 

Variable 


Column   Var.    Var . 
No .       type   name 

Line  1  of  group 

Group  repeated  ITIL  times 

1-5       R      RR1(K) 


description 


6-10 


R 


PS1(K) 


Random  roughness  after  Kth 
tillage  event  (0.0  -  4.0  cm) 

Soil  reflection  coefficient 
after  Kth  tillage  event  (0.0 
1.0). 


Line  2  of  group 

This  line  repeated  for  each  soil  horizon. 


1-5 

R 

BDl(K.KD) 

6-10 

R 

CK(K,KD) 

11-15 

R 

CB(K,KD) 

16-20 

R 

CAE(K,KD) 

21-25 

R 

THW(K.KD) 

26-30 


Bulk  density  (g/cmJ) . 

Saturated  hydraulic  conductivity 
(cm/day) . 

Campbell  B  coefficient. 

Campbell  air  entry  potential  (cm) 


THW(K,KD)    Water  content  at  permanent  wilt 
3    3 
(cm  /cm  ) . 

THS(K,KD)    Water  content  at  saturation 

3    3 

(cm  /cmJ) .   (K  =  subscript  for 

tillage  event,  KD  =  subscript 
for  soil  horizon. ) 
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Data  Group  19- 

Fertilizer 

Applications 


(One  line  must  be  present  even  if  blank;  insert  all 
fertilizer  application  lines  as  a  single  group.) 


Column   Var  • 
No.       type 


Var 


name 


Variable 
description 


Line  1  of  NF 


16 


All  units  in  lb/acre  of  species  shown  unless  otherwise 
labeled. 

1-5       R      FERT(l)      Depth  (cm)  of  a  uniform 

fertilizer  application. 
If  =  0.0,  a  surface  application  is 
indicated. 

6-10     R      FERT(2)      NHj  (0.0-400.0  lb/acre). 

11-15    R      FERT(3)     NO3  (0.0-1,600.0  lb/acre). 

0 

II 
Urea  (NH2-C-NH2) (0 . 0-800 . 0 
lb/acre) . 

Ca++  (0.0-20,000.0  lb/acre). 

SO4  (0.0-20,000.0  lb/acre). 

CO3  (0.0-20,000.0  lb/acre). 

Fraction   of  N15    in  NH^-N 
(N15/(N14+N15)).     (Molina   NT 
model   only) . 

41-45  R  FERT(9)  Fraction   of  N15    in  NO3-N. 

(Molina  NT  model   only) . 


16-20 

R 

FERT(4) 

21-25 

R 

FERT(5) 

26-30 

R 

FERT(6) 

31-35 

R 

FERT(7) 

36-40 

R 

FERT(8) 

l^NF   is   the   number   of    fertilizer   appli- 
cations.     The  maximum  number   allowed   in   the 
program   is    25. 
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Data  Group  20-- 
Organic  Residue 
Applications 


(Omit  if  NBYPAS  =  1;  see  Data  Group  2;  insert  all 
organic  matter  application  lines  as  a  single  group.) 


Column   Var . 
No .       type 

Line  of  NO17 


Var. 


name 


Variable 
description 


1-5 


6-10 


11-15 


16-20 


21-25 


26-30 


31-35 


36-40 


41-45 


46-50 


51-55 


R 


R 


R 


R 


R 


R 


R 


DEPTH 


ACN1 


SSAMT 


DELX(l) 


PR1 


AICODE 


AKLAB 


AKRES8 


APER 


AID 


Depth  (cm)  of  residue  application 
(equals  0  for  surface  applications) . 

Carbon: nitrogen  ratio  of  residue 
added  (for  example,  40:1  would 
be  entered  as  40.0). 

Organic  nitrogen  added  (lb/acre)  as 
oven  dry  weight  (0.0  -  100.0 
lb/acre) . 

Depth  (cm)  of  residue  on  the 
surface  (0-10cm) . 

Fresh  residue  reflection 
coefficient  (0.0-1.0) . 

Residue  identification  code  (see 
table  3) . 

Specific  rate  coefficient  for 
labile  fraction  (see  table  3) . 

Specific  rate  coefficient  for 
resistant  fraction  (see  table  3). 

Fraction  of  total  residue  that 
is  labile  (see  table  3) . 

Fraction  of  N15  in  residue  nitrogen 
application  (N15/(N14+N15) )   (Molina 
NT  model  only. ) 

Pure  carbon  switch. 

If  =  0.0,  normal  residue  N 

application. 
If  =  1.0,  pure  carbon  is  being 

added  in  SSAMT,  (cols.  11-15), 

(lb  C/acre) .   (Molina  NT  model 

only) . 


I'No  is  the  number  of  organic  residue  appli- 
cations. The  maximum  number  allowed  in  the 
program  is  5. 
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Table  3. 

Residue  identification  codes 


AICODE 


Values-1-  for 


AKLAB' 


AKRESIS- 


0  Variable 

1  Corn  stover 

2  Wheat  straw 

3  Cellulose 

4  Soybean  straw 

5  Variable 


User 

User  or 

supplied 

default 

0.317 

U  or  D 

.293 

U  or  D 

.098 

U  or  D 

.566 

U  or  D 

(4) 

U  or  D 

User  or 
default 

U  or  D 

U  or  D 

U  or  D 

U  or  D 

U  or  D 


Default  values  are  activated  by  leaving 

the  data  field  blank. 
2Default  AKLAB  =  0.044. 
3Default  AKRESIS  =  0.00071. 
4S  =  0.07  +  1.11  (1.0/ACNl)1-0/3'0. 


Data  Group  21 --Plant   (Omit  if  IVARR  =  2;  see  Data  Group  2.) 
Root  Distribution 

Column   Var . 


No .       type 
Line  1  of  1 


1-10 


R 


Var 


name 


KP2(1) 


11-20 

R 

KP2(2) 

21-30 

R 

KP2(3) 

31-40 

R 

KP2(4) 

41-50 

R 

KP2(5) 

51-60 

R 

KP2(6) 

Variable 
description 

Fraction  (decimal)  of  plant  roots 
from  0-  to  1-foot  depth. 

Fraction  from  1  foot  to  2  feet. 

Fraction  from  2  feet  to  3  feet. 

Fraction  from  3  feet  to  4  feet. 

Fraction  from  4  feet  to  5  feet. 

Fractions  >  5  feet. 
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Data  Group  22- -Plant  (Omit  if  CROP 
Uptake  of  Nitrogen 

Column  Var . 

No .  type 

Line  1  of  NPlf 


4;  see  Data  Group  2.) 


Var. 
name 


Variable 
description 


Data  Group  23 — 
Crop  Type  and 
Coefficients 


1-10     X      none        Space  for  identification. 

11-20    R      UPTK(K)      Plant  uptake  of  nitrogen 

(lb-N/acre  semimonthly) .  (K  = 
subscript  for  each  semimonthly 
period. ) 

(Present  only  if  CROP  =  4;  see  Data  Group  2.) 


Column 

No  . 


Var. 
type 


Line  1  of  1 
1-2      I 


3-10 

11-15 
16-20 

21-25 

26-33 
34-38 
39-43 


R 

R 
R 

R 


R 


R 


Var. 


name 


ICID 


CLASS 

PPOP 

DDTA 

DDMA 

ICR0P8 
AACOEF 
C0EF1 


Variable 
description 


Crop  Identification  Code  (l=corn, 
2=sorghum,  3=sweet  corn,  4=soybeans , 
5=spring  wheat,  6=winter  wheat,  7= 
oats,  8=barley,  9=rye ,  10=sunf lowers , 
ll=potatoes ,  12=alfalfa,  13=pasture 
grass,  14=sugar  beets,  15=cotton, 
16=peanuts ,  17=tomatoes ,  18=peas ,  19= 
carrots,  20=sugarcane) . 

Crop  maturity  class  for  corn  (in 
days,  for  example,  85  =  85 -day  corn). 

Plant  population  (No .  plants/ha) . 

Growing  degree  days  to  flowering 
(optional:   model  will  estimate  if 
blank  or  zero) . 

Growing  degree  days  to  maturity 
(optional:   model  will  estimate  if 
blank  or  zero) . 

Crop  variety  (for  identification 
only) . 

Optional  crop  transpiration 
coefficient.   Default=0 . 0085 . 

Optional  coefficient  for  photo- 
synthesis.  Default=65 . 0 . 


L8 


NP  is  24. 
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Column   Var.    Var . 
No.      type    name 

44-48    R      C0EF2 


49-53 


54-51 


59-63 


64-68 


69-73 


R 


R 


R 


C0EF3 


C0EF4 


C0EF5 


C0EF6 


C0EF7 


Variable 
description 

Optional  tops  coefficient  for 
respiration.   Default=0 . 74 . 

Optional  tops  coefficient  for 
production  of  leaves  and  stems. 
Defaults:   corn=3.58,  soybeans= 
4.10. 

Optional  tops  coefficient  for 
production  of  grain  or  other 
product.   Defaults:   corn=3.58, 
soybeans=3 . 00 . 

Optional  tops  coefficient  for 
total  drymatter  production. 
Defaults:   corn=0.0065,  soybeans^ 
0.0035. 

Optional  root  coefficient  for 
root  branching  (early  stages) . 

Optional  root  coefficient  for 
root  branching  (middle  stages) . 


Data  Group  24- -Soil 
Surface  and  Location 
Data 


Column   Var . 


Var 


No. 

tvpe 

name 

Line  1 

of  3 

1-10 

R 

RR 

11-20 

R 

PS 

21-30 

R 

SLOPE 

31-40 

R 

SLOPEL 

41-50 

R 

BDCOR 

51-60 


SL0PA8 


Variable 
description 


Random  roughness  (0.0  -  4.0  cm). 

Soil  reflection  coefficient 
(0.0  -  1.0) . 

Percent  slope  (0.0  -  50.0). 

Slope  length  (meters) . 

Bulk  density  of  stabilized  plow 
layer  (g/cnr ) . 

Slope  aspect  (N,  S,  E,  W,  NE , 
NW,  SE,  SW) . 
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Column   Var . 
No .      type 

Line  2  of  3 


1-10 


R 


Line  3  of  3 


Data  Group  25- - 

Input -Output  Options  No 

and  Run  Conditions 


21-25 

X 

26-30 

X 

31-35 

X 

36-40 

X 

41-45 

X 

46-50 

X 

51-55 

I 

Var. 
name 


WIND 


1-10 

R 

LAT 

11-20 

R 

ELEV 

Column 

Var. 

Var. 

No. 

type 

name 

Line  1 

of  2 

1-5 

X 

none 

6-10 

X 

none 

11-15 

X 

none 

16-20 

I 

ITENTH 

none 


none 


none 


none 


none 


none 


I  PORT 


Variable 
description 


Long-term  average  total  wind 
(km/day) .   (Note  that  daily  wind 
values  can  also  be  input;  set 
IDAYP8  =  1 . ) 


Latitude  of  site  location  (degrees) , 
Elevation  of  site  (meters). 

Variable 
description 

Not  used  by  program. 
Not  used  by  program. 
Not  used  by  program. 

Option  to  allow  printed  results 
to  be  output  at  0.1 -day  intervals 

If  ITENTH  =  1,  print  results  every 
10th  day. 

Not  used  by  program. 

Not  used  by  program. 

Not  used  by  program. 

Not  used  by  program. 

Not  used  by  program. 

Not  used  by  program. 

Option  to  select  the  special  form 
(Campbell  method)  of  subroutine 
PROP  to  compute  conductivities 
and  diffusivities  as  a  function 
of  moisture  content. 

If  IPOPT  =  1,  then  use  this  form. 
Campbell  coefficients  must  then 
be  included. 


6  6 


II 


Column   Var . 
No .       type 

Line  2  of  2 

1-5      I 

6-10     I 


11-15 


Var, 


name 


AA 


BB 


CC 


16-20 

X 

none 

21-25 

X 

none 

26-30 

R 

BBC 

31-35 

X 

none 

36-40 

X 

none 

41-45 

X 

none 

46-50 

I 

M 

51-55  X      none 

56-60  X      none 

Data  Group  26--Soil   Column  Var.    Var. 

Physical  Properties   No.  type   name 

Line  1  of  variable 

1-5  I       ICRACK 


6-10 


11-15 


I CRUST 


ICRCST 


Variable 
description 


If  AA  =  1,  print  the  initial  input 
data. 

If  BB  =  1,  print  computed  water 
contents  on  a  daily  basis. 

If  CC  =  1,  read  initial  water 
contents  from  Group  26. 

Not  used  by  program. 

Not  used  by  program. 

Volumetric  water  content  at  the 
bottom  boundary  expressed  as  a 
decimal  (dimensionless) .   (Note: 
This  can  be  less  than  or  equal  to 
the  saturated  value . ) 

Not  used  by  program. 

Not  used  by  program. 

Not  used  by  program. 

Minimum  number  of  time  intervals 
per  day  for  applying  the  finite 
difference  solution. 

Not  used  by  program. 

Not  used  by  program. 

Variable 
description 


If  =  0,  soil  is  not  a  cracking  soil. 

If  =  1,  soil  is  a  cracking  soil. 

If  =  0,  soil  surface  does  not  crust 

If  =  1,  soil  surface  does  crust. 

If  =  0,  soil  crust  does  not  crack. 

If  =  1,  soil  crust  does  crack. 
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Column 

Var. 

Var. 

No. 

tvpe 

name 

16-20 

R 

ARAN 

21-25 

26-30 
31-35 

36-40 

41-45 
46-50 


R 

X 
R 

R 

R 

R 


ARACA 


none 


WSEAL 


THSEAL 


THCR 


ARACR 


Variable 
description 

Maximum  rate  of  infiltration 
(cm/day)  for  normal  soil. 

Maximum  rate  of  infiltration 
(cm/day)  for  cracked  soil. 

Not  used  by  program. 

Precipitation  amount  (cm)  on  an 
unsealed  soil  above  which  sealing 
occurs . 

Soil  water  content  (9)    (cm  /cm^) 
on  an  unsealed  soil  above  which 
sealing  occurs. 

Soil  water  content  (8)    (cmJ/cmJ) 
below  which  cracking  occurs . 

Maximum  infiltration  rate  (cm/day) 
for  crusted  soil. 


51-55     R      ARACR8      Maximum  infiltration  rate  (cm/day) 

for  a  cracked  crust. 


Line  2  of  variable  (required  if  ICRUST  =  1) 
1-10      R      CRKSAT(l) 


19 


11-20 


R 


BEMPCR 


Saturated  hydraulic  conductivity 
(cm/day)  for  the  crusted  condition. 

Campbell  B  coefficient  for  the 
crusted  condition. 


21-30    R      AIREN7      Air-entry  potential  (cm)  for  the 

crusted  condition. 


31-37 


38-44 


45-51 


R 


R 


R 


CRTS 


CRTD 


CRTHI 


Saturated  water  content  (9) 

3    3 
(cmJ/cmJ)  for  the  crusted  condition. 

Water  content  (9)    at  permanent 

3    3 
wilt  (cm  /cnr)  for  the  crusted 

condition. 

3  3 

Initial  water   content    (9)    (cmJ/cm   ) 

for    the    crusted   conditions. 


l°Crusted   condition   is   assumed    to   occur   in 
the   uppermost   soil  horizon. 
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Column 

Var. 

Var. 

No. 

type 

name 

52-58 

R 

CRSAND 

59-65 

R 

CRSILT 

66-72 

R 

CRCLAY 

Variable 
description 

Percent  sand  for  the  crusted 
condition. 


Percent  silt  for  the  crusted 
condition. 

Percent  clay  for  the  crusted 
condition. 


Line  3  of  variable  (required  if  ICRCST  =   1) 
R      CRKSAT(2) 


1-10 


Saturated  hydraulic  conductivity 
(cm/day)  for  the  cracked  crust 
condition. 

11-20     R      BEMPC8       Campbell  B  coefficient  for  the 

cracked  crust  condition. 

21-30     R      AIREC8       Air-entry  potential  (cm)  for  the 

cracked  crust  condition. 

Line  4  of  variable  (always  present).    indicates  "normal" 
soil 

1-10  R 

11-20  R 

21-30  R 

31-37  R 

38-44  R 


KSAT(J.l)         Saturated  hydraulic   conductivity 
(cm/day)    for   normal    soil. 

BEMP(J)  Campbell   B   coefficient    for   normal 

soil . 

AIRENT(J)         Air-entry  potential    (cm)    for   normal 
soil . 

TS(J)  Saturated  water   content    (0) 

(cm-ycm   )    for   normal   soil. 

TD(J)  Water   content    (0)    at   permanent 

wilt    (cm  -ycmy    for   normal    soil. 


20a  paired   sequence   of    lines    4   and    5   must   be 
repeated    for  each  horizon. 

21"Normal"    soil   defined   as   not    cracked   or 
crusted. 
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Column   Var 


Var. 


No .       type    name 
45-51     R      THI(J) 


52-58 

R 

SAND(J) 

59-65 

R 

SILT(J) 

66-72 

R 

CLAY(J) 

73-79 

R 

TSM(J) 

Variable 
description 

Initial  water  content  (8)    (at 
start  of  run)  (cm-ycm  )  f°r 
normal  soil. 

Percent  sand  for  normal  soil. 

Percent  silt  for  normal  soil. 

Percent  clay  for  normal  soil. 

Leaching  water  content  (cm  /cm  ) 
for  Jth  horizon.  (J  =  subscript 
for  Jth  horizon.) 

22 


Line    5    of  variable    (required    if   ICRACK  =   1) 

1-10  R  KSAT(J,2)         Saturated  hydraulic    conductivity 

(cm/day)    for   cracked   condition. 

11-20  R  BEMPCA(J)         Campbell    B   coefficient    for 

cracked  condition. 

21-30  R  AIREN8(J)         Air-entry  potential    (cm)    for 

cracked  condition.       (J   = 
subscript   for  Jth  horizon.) 


"A  paired   sequence   of   lines    4   and   5  must   be 
repeated   for  each  horizon. 
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WATER  APPLICATION  FILE 


This  input  file  can  be  generated  with  the  help  of  the 
interactive  operating  system  (see  "Interactive  Generation  of 
Input  Files").   Otherwise,  a  suitable  text  editor  or  keypunch 
can  be  used  to  follow  the  formats  as  shown. 


Column   Var . 
No.       type 

Line  1  of  set 

11-13     I 

Line  2  of  set 


ni-12 

13-14 
15 


1 
I 

A 

R 


16-21 
**22-26  X 

27-28     I 


Var. 
name 


APPS 


MON(I) 
DATE ( I ) 
ADENT(I) 

AMT(I) 

none 

M0N(I+1) 


Variable 
description 


Number  of  water  applications 
(max  =  99) . 


Month. 

Day. 

Type  of  application,  R  =  rain, 
I  =  irrigation. 

Depth  of  application  (inches). 

Not  used  by  program. 

Month. 


Repeat  from  *  to  **  three  times  per  line;  last  entry  will  be  for 
columns  64-69  and  will  contain  depth  of  application.  Additional 
lines  can  be  added  if  necessary.  Maximum  number  of  applications 
is  99. 
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PAN  EVAPORATION,  WIND,  AND  SOLAR 
RADIATION  FILE 


This  particular  input  file  has  two  possible  structures  depending 
on  whether  the  input  data  are  daily  or  semimonthly.   Both 
formats  can  be  generated  with  the  assistance  of  the  interactive 
operating  system  (see  "Interactive  Generation  of  Input  Files"). 

For  semimonthly  inputs : 


Column   Var . 

Var. 

Variable 

No .       type 

name 

description 

Line  1  of  2 

11-15    R      Ul  Pan  evaporation  (inches)  for  Jan. 

1-15. 

16-20    R      Ul  Pan  evaporation  (inches)  for  Jan. 

16 -end  of  month. 

21-25    R      Ul  Pan  evaporation  (inches)  for  Feb. 

1-15. 


Repeat  in  semimonthly  increments  to: 

66-70     R      Ul  Pan  evaporation  (inches)  for  June 

16 -end  of  month. 

Line  2  of  2 

11-15    R      Ul  Pan  evaporation  (inches)  for  July 

1-15. 


Repeat  in  semimonthly  increments  to : 

66-70     R      Ul  Pan  evaporation  (inches)  for  Dec. 

16 -end  of  month. 

Note:   Wind  and  solar  radiation  data  are  not  included  if 
semimonthly  inputs  are  selected.   This  option  may  be  added  in 
the  future . 


7  2 


For  daily  inputs: 

Column   Var .    Var . 
No .       type    name 

Line  1  of  N1 

*l-5     R      PAN 


6-10     R 


**11-15   R 


16-20    R 


WIND 


SOLAR 


PAN 


Variable 
description 


Pan  evaporation  (inches)  for  first 
day  of  week. 

Wind  run  (km)  for  first  day  of 
week. 

Solar  radiation  (langleys)  for 
first  day  of  week. 

Pan  evaporation  (inches)  for  second 
day  of  week. 


Repeat  from  *  to  **  for  each  day  of  week;  last  entry  will  be  in 
columns  101-105. 

A  line  is  needed  for  each  week  of  data  included.   First  day 
should  correspond  to  starting  day  of  run. 


_N  is  the  number  of  weeks . 
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APPENDIX  A- -GENERAL  SUMMARY  OF  OUTPUT  FILES  AND  SPECIFIC  FILE 
STRUCTURES 


FILE 
Output 


CONTENTS 


Summary  of  general  control 
parameters  and  model  input 
data. 

Predicted  soil  profile 
data  for  major  cations  and 
nitrogenous  species 
at  specified  times. 
Mass  balance  status  for 
nitrogenous  species  at 
specified  times . 
Miscellaneous  monitoring 
information. 


LU  2 


LU  4 


LU  6 


LU  7 


Soil  leachate  data  for  flow 
volumes  and  soluble 
constitute  masses  becomes 
input  to  aquif ier  submodels . 

Daily  output  from  crop 
growth  model: 

1.  Total  dry  matter 
(tops)  (kg/ha). 

2.  Leaf  area  (dm  /plant) . 

3.  Corn  grain  (kg/ha). 

4.  N  uptake  demand  (kg/ha). 

Daily  output  from  root 
growth  model : 

1.  Total  net  root  length 

o 
produced  (cm  root/cm 

soil) . 

2.  Total  root  length 
produced  (cm  root/cm 
soil) . 

3 .  Total  root  length 
sloughed  (cm  root/cm^ 
soil) . 

4.  Net  root  mass  produced. 

5.  Root  status  (0-1  range). 

Daily  summary  of  root 
distribution  in  each  soil 
segment  (for  example,  15  cm) 
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LU  11 


LU  12 


Daily  summary  from 
unsaturated  flow  model. 

1.  Basic  monitoring 
parameters  including 
cumulative  leachate  (cm) 
and  evapotranspiration  (cm) 

2 .  Soil  water  contents 
(cm  /cm  )  at  each  node. 

3.  Soil  matrix  potential  (cm 
HoO)  at  each  node. 

Output  summary  from  nitrogen 
transformation  submodel 
(Molina) . 


LU  13 


Root  densities,  soil  water 

3    3 
contents  (9)    (cm  /cm  ) ;  soil 

matrix  potential  (VO  (cm 

H2O) ;  soil  temperatures  (°C) 

as  a  function  of  time  and 

depth  (cm)  (suitable  for  3-D 

plots) . 


LU  14 


Daily  summary  of  surface 
conditions . 


LU  15 


Summary  of  N15  output  from 
Molina  submodel. 


LU  17 


Summary  of  daily  evapo- 
transpiration. 


Unsaturated  Flow 
Input  (LU  1) 


Record   Type     Format     Variable     Contents 

1      Binary   NA         II  Dummy. 

1      Binary   NA         12  Dummy. 

1      Binary   NA         13  Dummy. 

1      Binary   NA         13  Dummy. 

1       Binary    NA         CMH201(1)     Volume  of  water  in 

segment  1  (cm/cm^). 
1       Binary    NA         M01S1N(1)     Water  flux  into 

segment  1  (era/cm^) . 
1      Binary   NA         M01S08(1)     Water  flux  out  of 

segment  1  (cm/cm  ) . 
1      Binary   NA         TEN1(1)       Soil  suction  for 

segment  1  (+cm  H2O) 
1       Binary    NA         U(l)  Transpiration  for 

segment  1  (cm/cm^) . 

Entire  record  sequence  is  repeated  for  each  soil  segment 
starting  at  the  soil  surface.   NA  =  not  applicable. 
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Soil  Leachate  Data 
(Input  to  Aquifer 
Submodels)  (LU  2) 


Record    Type 


Coded 
Coded 
Coded 

Coded 

Coded 


Format 

(IX,  15, 

15, 
F10. 2, 

F10. 2, 

F10. 2, 


Coded     F10. 2 


Coded     F10. 2 


Coded    F10. 2 


Coded    F10. 2 


Coded    F10. 2 


Variable 

YEAR 
I  DAY 
DELN 

CPPM(l) 

CPPM(2) 

CPPM(3) 

CPPM(4) 

CPPM(5) 


Coded     F10. 2,      CPPM(6) 


CPPM(7) 


Coded     F10. 2,      CPPM(8) 


CPPM(9) 


Coded     F10. 2)      CPPM(IO) 


Contents 

Year  number . 
Julian  day. 
Leachate  for  this 

time  step  (cm) . 
Cone.  NO3-N  (mg/L) 

leached. 
Cone.  NH^-N  (mg/L) 

leached  for  this 

time  step. 
Cone.  UREA-N 

(mg/L)  leached 

for  this  time  step 
Cone.  Ca++  (mg/L) 

leached  for  this 

time  step. 
Cone.  Na+  (mg/L) 

leached  for  this 

time  step. 
Cone.  Mg++  (mg/L) 

leached  for  this 

time  step. 
Cone.  HCO3  (mg/L) 

leached  for  this 

time  step. 
Cone.  CI"  (mg/L) 

leached  for  this 

time  step. 
Cone.  C0|  (mg/L) 

leached  for  this 

time  step. 
Cone.  SO4  (mg/L) 

leached  for  this 

time  step. 


Daily  Output  From 
Crop  Growth  Model 
(LU  4) 


Record 

1 
1 

1 
1 
1 


Type 

Coded 
Coded 

Coded 
Coded 
Coded 


Format 


(7X 

18X 


13, 
E9.4 


Variable 

I 
PDRY 


15X.E9.3,  ALAD 
19X.E9.4,  PEAR 
21X.E8.3)   TNUPT8 


Contents 

Julian  day  number. 
Total  dry  matter 

(tops)  (kg/ha). 
Leaf  area  (dm^/plant) 
Corn  grain  (kg/ha) . 
N  uptake  demand 

(kg/ha) . 


Repeat  record  sequence  for  each  Julian  day. 
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Daily  Output  From 
Root  Growth  Model 
(LU  6) 


Record 

1 
1 


Type 

Coded 
Coded 

Coded 


Format 


Variable 


(7X.I5,     IDAY 
11X.F10.3,  RRTOT 

10X.F12.3,  TDM 


Coded     12X.F12.3,  TDEAD 

Coded    12X.F12.3,  ROOTM 
Coded    12X.F6.3)   RATIO 


Contents 

Julian  day  number. 
Total  net  root  Length 

(cm  root/cm^  soil) 
Total  root  length 

produced  (cm  root/ 

cm^-  soil)  . 
Total  root  length 

sloughed  (cm  root/ 

cm^-  soil)  . 
Net  root  mass  produced 

(kg/ha) . 
Root  status  (0-1 
range) . 


Repeat  record  sequence  for  each  Julian  day. 


Daily  Summary  of 
Root  Distribution 
(LU  7) 


Record 

Type 

Format 

Variable 

Contents 

1 

Coded 

(9X.I5 

IDAY 

Julian  day  number. 

1 

Coded 

7X.F5.0, 

DELX(2) 

DELX  for  segment  2 . 

1 

Coded 

12X.F7.3, 

KP2(1) 

Fraction  of  roots  in 
segment  2 . 

1 

Coded 

F7.3, 

DELX(3) 

DELX  for  segment  3 . 

1 

Coded 

F7.3, 

KP2(2) 

Fraction  of  roots  in 
segment  3 . 

Repeat  to : 

1      Coded 

F7.3, 

DELX(ll) 

DELX  for  segment  11. 

1      Coded 

F7.3) 

KP2(10) 

Fraction  of  roots  in 
segment  10. 

Repeat  record  sequence  for  each  Julian  day, 
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Output  Summary  From  Record 
Nitrogen  Trans- 
formation Submodel  1 
(Molina)  (LU  12)  1 


Type 

Format 

Variable 

Contents 

Coded 

(IX, 15 

I  DAY 

Julian  day  number. 

Coded 

F7.1, 

TAMNI 

Total  NHj-N  +  NO3-N 
in  soil  (kg/ha) . 

Coded 

F7.1, 

DPLT 

Daily  plant  uptake  of 
N  (kg/ha). 

Coded 

F7.1, 

DOUT 

Daily  N  leached 
(kg/ha) . 

Coded 

F7.1, 

DDENIT 

Daily  N  denitrified 
(kg/ha) . 

Coded 

F7.1, 

SUMR11/10 

Daily  N  fixed  (kg/ha) 
--N  on  surface 
(kg/ha) . 

Coded 

F7.1, 

AN03(1)/10 

NO3-N  on  surface 
(kg/ha) . 

Coded 

F7.1) 

ANH3(1)/10 

NHj-N  on  surface 
(kg/ha) . 

Repeat  record  sequence  for  each  Julian  day. 


Summary  of  Root  and 
Soil  Data  for  3-D 
Plotting  (LU  13)1 


Record 

1 
1 
1 

1 
1 


Type 

Coded 
Coded 
Coded 

Coded 
Coded 

Coded 

Coded 

Coded 
Coded 


Format 

(IX, 15, 

15, 
F8.3, 

F8.3, 
F8.3, 

F10. 3, 

F10. 3, 

F10. 3) 
(F10.3, 


Variable 

I  DAY 

JDEP1 

RT0T10/DELX 

TEM 
CMH201/DELX 

TEN1 

BD 

AN03/CMH201 

ANA/CMH201 

ACA/CMH201 

AMG/CMH201 

HC03/CMH201 

CL/CMH201 

S04/CMH201 


Contents 

Julian  day  number. 

Depth  (cm) . 

Root  length  (cm  root/ 

cm  soil) . 
Soil  temperature  (°C) 

Soil  water  content 

3    3 

(cmJ/cmJ) . 

Soil  suction  (cm 

water) . 
Soil  bulk  density 

(g/cm3) . 

NO3-N  cone.  (mg/L) . 

Na  cone.  (mg/L). 
,++ 


cone 


Mg    cone 


CI"  cone. 
SO^  cone . 


(mg/L) 
(mg/L) 
(mg/L) 

(mg/L) . 

(mg/L) . 


^Entire  file  is  shown;  if  option  1  is 
selected  for  IWRI3,  salts  are  not 
included. 


Repeat  this  sequence  for  each  depth  (segment)  starting  near  the 
soil  surface. 


Repeat  the  entire  record  for  each  Julian  day- 


Daily  Summary  of 
Surface  Conditions 
(LU  14) 


Record   Type 


Coded 
Coded 
Coded 

Coded 

Coded 


Format     Variable      Contents 

(15,        IDAY         Julian  day  number. 
F10 . 3 ,      RR  Random  roughness  (cm). 

F10. 3,      BD(2)        Bulk  density  segment  2 

(g/cm3). 
F10. 3,      SURRES        Surface  residue 

(kg/ha) . 
F10. 3)      PCOVER       Percent  cover. 


Repeat  record  sequence  for  each  Julian  day. 


Summary  of  N15 
Output  From  Molina 
Submodel  (LU  15) 


Record 

Type 

Format 

Variabl 

e 

Contents 

1 

Coded 

(5X.49H, 

Holleri 

th 

Heading. 

1 

Coded 

15) 

IDAY 

Julian  day  number 

2 

Coded 

(5X.49H, 

Holleri 

th 

Heading. 

2 

Coded 

5X.21H) 

Holleri 

th 

Heading. 

3 

Coded 

(3X.I5, 

I 

Segment  number. 

3 

Coded 

F13.3, 

PER8 

TNHt-N  (%  N15) . 

3 

Coded 

F18.3, 

PER8S 

SNHJ-N  (%  N15) . 

3 

Coded 

F14.3, 

PERU 

NO3-N  (%  N15) . 

3 

Coded 

F14.3, 

PERMM 

Microbial  mass 
(%  N15) . 

3 

Coded 

F14.3) 

PERNH 

Nutrient  humus 
(%  N15) . 

Repeat  record  sequence  for  each  Julian  day  and  depth. 
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Summary  of  Daily 
Evapo transpiration 
Output  (LU  17) 


Record 

Type 

Format 

Variable 

Contents 

L 

Coded 

(IX,  12, 

II 

Julian  day  number. 

1 

Coded 

E12.4, 

DEVE 

Surface  evaporation 
(cm/mo) . 

1 

Coded 

E12.4, 

DTRAN 

Transpiration  (cm/mo) 

1 

Coded 

E12.4, 

DAYET 

Evapo transpiration 
(cm/mo) . 

1 

Coded 

E12.4, 

DEVAP 

Potential  surface 
evaporation  (cm/mo) 

1 

Coded 

E12.4, 

POTRAN 

Potential  transpira- 
tion (cm/mo) . 

1 

Coded 

E12.4, 

POET 

Potential  evapo - 
transpiration 
(cm/mo) . 

1 

Coded 

E12.4, 

1-TR 

Fraction  of  radiation 
reaching  ground. 

1 

Coded 

E12.4) 

UDAY 

Pan  evaporation 
(cm/mo) . 

Repeat  record  sequence  for  each  Julian  day. 
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APPENDIX  B- -SAMPLE  RUNS 


Sample  Interactive 
Session  To  Create 
Input  Files  for  Soil 
Temperature  Submodel 


GREETINGS  --  PROCEDURE  TO  RUN  NTRM  SOIL  TEMPERATURE 
SUBMODEL 


ENTER  CONTROL  CODE: 

1=  CREATE  MAX-MIN  AIR  TEMPERATURE  FILE 
2=  CREATE  GENERAL  SOIL  TEMPERATURE  INPUT  FILE 
3=  CREATE  BOTH  FILES 
3 

ENTER  NAME  OF  NEW  MAX-MIN  AIR  TEMP  FILE: 
NEWTEMP 


ROUTINE  TO  GENERATE  MAX-MIN  AIR  TEMPERATURE  FILE 

PLEASE  ENTER  YOUR  JULIAN  STARTING  DAY 
121 


ENTER  MAX  AND  MIN  AIR  TEMPERATURES  FOR  EACH  DAY. 
FREE  FORMAT. 

ENTER  -999,-999  WHEN  FINISHED.   UNITS  ARE  DEGREES 
EXAMPLE:     56,45 


USE 


71, 

40 

70, 

57 

65, 

54 

59, 

52 

60, 

47 

58, 

44 

60, 

39 

59 

47 

51 

39 

54 

35 

60 

35 

60 

51 

67 

41 

73 

41 

77 

43 

74 

46 

61 

54 

71 

45 

73 

39 

78 

44 

81 

47 

77 

59 

70 

60 

62 

,54 

59 

,52 

64 

,60 

71 

,54 

69 

,56 

77,59 

70,49  [ 

76,42  L 

80,55 

77.54  r 
74,53  I. 
78,49 

83,60  s 

83,53 

87,57  L 

70,56 

72.55  f 
-999,-999 


MAX  -  MIN  AIR  TEMPERATURE  FILE  CREATED 


ENTER  NAME  OF  NEW  GENERAL  INPUT  FILE  FOR  SOIL  TEMP: 
NESOTEM 


ROUTINE  TO  SET  UP  GENERAL  SOIL  TEMPERATURE  INPUT  FILE. 


ENTER  YOUR  JULIAN  STARTING  DAY  AND  NUMBER  OF  DAYS  IN  RUN. 
USE  FREE  FORMAT. 


121,40 


ENTER  INITIAL  (AT  TIME  ZERO  HOUR)  SOIL  TEMPERATURE 
AT  EACH  OF  THE  INDICATED  NODE  DEPTHS  (CM). 
USE  FREE  FORMAT.   UNITS  ARE  DEGREES  C. 

NODAL  DEPTH  (CM) 

.1     .2     .5    1.0    1.5    2.0    5.0   10.0   15.0   20.0 
15,15,15,15,15,15,15,15,15,15 

NODAL  DEPTH  (CM) 

25.    30.    35.    40.    45.    50.    55.    60.    65.    70. 
15,15,15,15,15,15,15,15,15,15 

NODAL  DEPTH  (CM) 

75.    80.    85.    90.    95.   100.   110.   120.   135.   150. 
16,16,17,17,17,18,18,18,19,19 
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NODAL  DEPTH  (CM) 

200.   250.   300.   350.   400.   450.   500.   550 
20,20,20,20,20,20,20,20 


ENTER  YOUR  MANAGEMENT  CODE: 

1  =  BARE  SOIL 

3  -  50  PCT.  OAT  RESIDUE  COVER 
5  -  75  PCT.  OAT  RESIDUE  COVER 
7  =  99  PCT.  OAT  RESIDUE  COVER 
9  =  CORN  +  50  PCT.  OAT  COVER 
11  =  CORN  +  NO  RESIDUE  COVER 
11 

GENERAL  SOIL  TEMPERATURE  INPUT  FILE  CREATED 


Contents  of  Example   121  71.0  40.0 
Max-Min  Air  122  70.0  57.0 

Temperature  File      123  65.0  54.0 

124  59.0  52.0 

125  60.0  47.0 

126  58.0  44.0 

127  60.0  39.0 

128  59.0  47.0 

129  51.0  39.0 

130  54.0  35.0 

131  60.0  35.0 

132  60.0  51.0 

133  67.0  41.0 

134  73.0  41.0 

135  77.0  43.0 

136  74.0  46.0 

137  61.0  54.0 

138  71.0  45.0 

139  73.0  39.0 

140  78.0  44.0 

141  81.0  47.0 

142  77.0  59.0 

143  70.0  60.0 

144  62.0  54.0 
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145  59.0  52.0 

146  64.0  60.0 

147  71.0  54.0 

148  69.0  56.0 

149  77.0  59.0 

150  70.0  49.0 

151  76.0  42.0 

152  80.0  55.0 

153  77.0  54.0 

154  74.0  53.0 

155  78.0  49.0 

156  83.0  60.0 

157  83.0  53.0 

158  87.0  57.0 

159  70.0  56.0 

160  72.0  55.0 


Contents  of  Example   121 
General  Input  File     38 

for  Soil  Temperature    .1   15.0   .2   15.0   .5   15.01.0   15. 01. 5   15.02.0   15.0 
Submodel  5.0   15.010.0   15.015.0   15.020.0   15.025.0   15.030.0   15.035.0 

15.040.0   15.0 

45.0  15.050.0   15.055.0   15.060.0   15.065.0   15.070.0   15.075.0 
16.080.0   16.0 

85.   17.0  90.   17.0  95.   17.0100.   18.0110.   18.0120.   18.0135. 
19.0150.   19.0 

200.   20.0250.   20.0300.   20.0350.   20.0400.   20.0450.   20.0500. 
20.0550.   20.0 
40 
1 
11   1 
13.9   13.9   13.9   13.9   13.9   13.9 

13.9   13.9   13.9   18.4   18.4   18.4   22.1   22.1   22.1   27.7 
27.7   27.7   27.5 

27.5   27.5   24.4   24.4   24.4   26.8   26.8   26.8   25.2   22.6 

22.1  22.1   24.4 

24.4   24.4   24.4   24.4   24.4   24.4 
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Sample  Interactive 
Session  To  Create 
Common  Input  File 
for  NTRM  Central 
Model 


GENERATION  SYSTEM  FOR  COMMON  INPUT  FILE (DECK) 

ENTER  FILE  NAME  FOR  NEW  COMMON  INPUT  FILE  TO  BE  CREATED: 

NEWCOM 

THE  FOLLOWING  ARE  THE  INPUT  OPTIONS  PRESENTLY  AVAILABLE 

IN  THIS  OPERATING  PROCEDURE 


1.  DEFAULT  OPTION  1 

2.  DEFAULT  OPTION  2 


MINIMAL  USER  INPUT  REQUIRED 
ALLOWS  SELECTION  OF  SUBMODELS 
UTILIZED,  MORE  USER  INPUT  MAY 
BE  REQUIRED 


ENTER  OPTION  SELECTION 


INPUT  TITLE  DESIRED  ON  PRINTED  OUTPUT 
(1  TO  80  CHARACTERS) 


ROSEMOUNT  1981  --  TILL  CASE 


121 
160 


THE  FOLLOWING  INPUT  SETS  TIME  FRAMEWORK  FOR  THE 

SIMULATION  RUN 

INPUT  JULIAN  DATES  FOR  THE  FOLLOWING 

STARTING  DATE  OF  THE  RUN  - 

TERMINATION  DATE  OF  THE  RUN  = 


150 


ENTER  JULIAN  DATE  FOR  CROP  EMERGENCE 

THE  FOLLOWING  INPUT  SETS  SOIL   DATA 
INPUT  NUMBER  OF  SOIL  HORIZONS  (1  TO  10) 


15 


INPUT  DEPTH  (CM)  FROM  SURFACE  TO  BOTTOM  OF 
HORIZON  1 


ENTER  NUMBER  OF  SEGMENTS  FOR  HORIZON  1 


ENTER  NUMBER  OF  NODAL  SUBDIVISIONS  IN  EACH 
SEGMENT  (FOR  HORIZON  1) 


INPUT  DEPTH  (CM)  FROM  SURFACE  TO  BOTTOM  OF 
HORIZON  2 


45 


85 


105 


150 


.001 

3.4 

0 

.95 
.34 
L0 


ENTER  NUMBER  OF  SEGMENTS  FOR  HORIZON  2 


ENTER  NUMBER  OF  NODAL  SUBDIVISIONS  IN  EACH 
SEGMENT  (FOR  HORIZON  2) 


INPUT  DEPTH  (CM)  FROM  SURFACE  TO  BOTTOM  OF 
HORIZON  3 


ENTER  NUMBER  OF  SEGMENTS  FOR  HORIZON  3 


ENTER  NUMBER  OF  NODAL  SUBDIVISIONS  IN  EACH 
SEGMENT  (FOR  HORIZON  3) 


INPUT  DEPTH  (CM)  FROM  SURFACE  TO  BOTTOM  OF 
HORIZON  4 


ENTER  NUMBER  OF  SEGMENTS  FOR  HORIZON  4 


ENTER  NUMBER  OF  NODAL  SUBDIVISIONS  IN  EACH 
SEGMENT  (FOR  HORIZON  4) 


INPUT  THE  FOLLOWING  SOIL  CHEMICAL  ANALYSIS  FOR 
CHEMICAL  HORIZON     1 

ENTER  NH4+  (SOLUBLE  FORM)  (MEQ/L) 

(WARNING! !   NEVER  USE  TOTAL  NH4+  HERE) 

ENTER  N03-   (MEQ/L) 

ENTER  UREA   (MEQ/L) 

ENTER  BULK  DENSITY  (G/CM3  OF  SOIL) 

ENTER  GM  WATER/GM  SOIL  IN  SOIL  EXTRACT 

ENTER  MICROBIAL  BIOMASS  (UG/G  SOIL)  AS  CARBON 
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500 


.01 
2.4 

0 

.95 
.34 
5 
196 


.01 
.60 

0 

.95 
.34 

1 
24 


.01 


ENTER  NUTRIENT  HUMUS  (UG/G  SOIL)  AS  CARBON 

INPUT  THE  FOLLOWING  SOIL  CHEMICAL  ANALYSIS  FOR 
CHEMICAL  HORIZON     2 

ENTER  NH4+  (SOLUBLE  FORM)  (MEQ/L) 

(WARNING! !   NEVER  USE  TOTAL  NH4+  HERE) 

ENTER  N03-   (MEQ/L) 

ENTER  UREA   (MEQ/L) 

ENTER  BULK  DENSITY  (G/CM3  OF  SOIL) 

ENTER  GM  WATER/GM  SOIL  IN  SOIL  EXTRACT 

ENTER  MICROBIAL  BIOMASS  (UG/G  SOIL)  AS  CARBON 

ENTER  NUTRIENT  HUMUS  (UG/G  SOIL)  AS  CARBON 

INPUT  THE  FOLLOWING  SOIL  CHEMICAL  ANALYSIS  FOR 
CHEMICAL  HORIZON     3 

ENTER  NH4+  (SOLUBLE  FORM)  (MEQ/L) 

(WARNING!!   NEVER  USE  TOTAL  NH4+  HERE) 

ENTER  N03-   (MEQ/L) 

ENTER  UREA   (MEQ/L) 

ENTER  BULK  DENSITY  (G/CM3  OF  SOIL) 

ENTER  GM  WATER/GM  SOIL  IN  SOIL  EXTRACT 

ENTER  MICROBIAL  BIOMASS  (UG/G  SOIL)  AS  CARBON 

ENTER  NUTRIENT  HUMUS  (UG/G  SOIL)  AS  CARBON 

INPUT  THE  FOLLOWING  SOIL  CHEMICAL  ANALYSIS  FOR 
CHEMICAL  HORIZON     4 

ENTER  NH4+  (SOLUBLE  FORM)  (MEQ/L) 

(WARNING! !   NEVER  USE  TOTAL  NH4+  HERE) 
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.60 

0 

.95 

.35 

1 

2.3 

20 


20 


11.3 


15.3 


.30 


08 


22 


ENTER  N03-   (MEQ/L) 

ENTER  UREA   (MEQ/L) 

ENTER  BULK  DENSITY  (G/CM3  OF  SOIL) 
ENTER  GM  WATER/GM  SOIL  IN  SOIL  EXTRACT 
ENTER  MICROBIAL  BIOMASS  (UG/G  SOIL)  AS  CARBON 
ENTER  NUTRIENT  HUMUS  (UG/G  SOIL)  AS  CARBON 

ENTER  MAXIMUM  SURFACE  INFILTRATION  RATE  (CM/DAY) 


THE  FOLLOWING  WILL  SET  UP  SOIL  INFORMATION, 
AND  THE  UNSATURATED  FLOW  PROPERTIES  FOR  EACH 
SOIL  HORIZON 

ENTER  SATURATED  HYDRAULIC  CONDUCTIVITY  (CM/DAY) 
FOR  HORIZON  1 


ENTER  CAMPBELL  B  COEFFICIENT  FOR  HORIZON  1 


ENTER  AIR- ENTRY  POTENTIAL  (CM)  FOR  HORIZON  1 


ENTER  SATURATED  WATER  CONTENT  (CM3/CM3)  FOR  HORIZON  1 


ENTER  SOIL  WATER  CONTENT  (CM3/CM3)  AT  PERMANENT 


WILT  (AT  15  BARS)  FOR  HORIZON  1 


ENTER  THE  INITIAL  SOIL  WATER  CONTENT  (CM3/CM3)  FOR 
HORIZON  1 
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20 


20 


.22 


20 


11.3 


15.3 


30 


.08 


.22 


20 


20 


.22 


ENTER  PERCENT  SAND  FOR  HORIZON  1 


ENTER  PERCENT  SILT  FOR  HORIZON  1 


ENTER  LEACHING  WATER  CONTENT  (CM3/CM3)  FOR  HORIZON  1 


ENTER  SATURATED  HYDRAULIC  CONDUCTIVITY  (CM/DAY) 
FOR  HORIZON  2 


ENTER  CAMPBELL  B  COEFFICIENT  FOR  HORIZON  2 


ENTER  AIR- ENTRY  POTENTIAL  (CM)  FOR  HORIZON  2 


ENTER  SATURATED  WATER  CONTENT  (CM3/CM3)  FOR  HORIZON  2 


ENTER  SOIL  WATER  CONTENT  (CM3/CM3)  AT  PERMANENT 
WILT  (AT  15  BARS)  FOR  HORIZON  2 


ENTER  THE  INITIAL  SOIL  WATER  CONTENT  (CM3/CM3)  FOR 
HORIZON  2 


ENTER  PERCENT  SAND  FOR  HORIZON  2 


ENTER  PERCENT  SILT  FOR  HORIZON  2 


ENTER  LEACHING  WATER  CONTENT  (CM3/CM3)  FOR  HORIZON  2 
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20 


11.3 


15.3 


30 


08 


22 


20 


20 


.22 


20 


11.3 


15.3 


ENTER  SATURATED  HYDRAULIC  CONDUCTIVITY  (CM/DAY) 
FOR  HORIZON  3 


ENTER  CAMPBELL  B  COEFFICIENT  FOR  HORIZON  3 


ENTER  AIR- ENTRY  POTENTIAL  (CM)  FOR  HORIZON  3 


ENTER  SATURATED  WATER  CONTENT  (CM3/CM3)  FOR  HORIZON  3 


ENTER  SOIL  WATER  CONTENT  (CM3/CM3)  AT  PERMANENT 
WILT  (AT  15  BARS)  FOR  HORIZON  3 


ENTER  THE  INITIAL  SOIL  WATER  CONTENT  (CM3/CM3)  FOR 
HORIZON  3 


ENTER  PERCENT  SAND  FOR  HORIZON  3 


ENTER  PERCENT  SILT  FOR  HORIZON  3 


ENTER  LEACHING  WATER  CONTENT  (CM3/CM3)  FOR  HORIZON  3 


ENTER  SATURATED  HYDRAULIC  CONDUCTIVITY  (CM/DAY) 
FOR  HORIZON  4 


ENTER  CAMPBELL  B  COEFFICIENT  FOR  HORIZON  4 


ENTER  AIR-ENTRY  POTENTIAL  (CM)  FOR  HORIZON  4 
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30 


.08 


.22 


20 


20 


22 


.22 


1.0 


5 


0.0 


0 


ENTER  SATURATED  WATER  CONTENT  (CM3/CM3)  FOR  HORIZON  4 


ENTER  SOIL  WATER  CONTENT  (CM3/CM3)  AT  PERMANENT 


WILT  (AT  15  BARS)  FOR  HORIZON  4 


ENTER  THE  INITIAL  SOIL  WATER  CONTENT  (CM3/CM3)  FOR 
HORIZON  4 


ENTER  PERCENT  SAND  FOR  HORIZON  4 


ENTER  PERCENT  SILT  FOR  HORIZON  4 


ENTER  LEACHING  WATER  CONTENT  (CM3/CM3)  FOR  HORIZON  4 


ENTER  VOLUMETRIC  WATER  CONTENT  AT  THE  BOTTOM  NODE 
EXPRESSED  AS  A  DECIMAL  (DIMENSIONLESS) 


THE  FOLLOWING   GUIDES  THE  INPUT  FOR 
THE  EVAPORATION  SUBMODEL 


ENTER  RANDOM  ROUGHNESS  (0.0  -  4.0  CM) 


ENTER  SOIL  REFLECTION  COEFFICIENT  (0.0  -  1.0) 


ENTER  PERCENT  SLOPE  (0.0  -  50.0) 


ENTER  SLOPE  LENGTH  (METERS) 
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1.2 


200 


41 


300 


0.0 


0.1 


0.0 


0.0 


0.0 


ENTER  STABILIZED  BULK  DENSITY  OF  SURFACE  HORIZON 
(0  <  VALUE  <  2.6  GM/CM3) 


ENTER  LONG-TERM  AVERAGE  WIND  TRAVEL  (KM/DAY) 


ENTER  LATITUDE  OF  LOCATION  (DEGREES) 


ENTER  ELEVATION  OF  LOCATION  (METERS) 


THE  FOLLOWING  WILL  GUIDE  IN  THE  INPUT 
OF  RAIN  AND  IRRIGATION  WATER  ANALYSIS 


DO  YOU  WISH  TO  SET  RAIN  WATER  ANALYSIS 
TO  DEFAULT  VALUE  OF  0.0.   ENTER  1  FOR  YES 
0  FOR  NO 


IS  IRRIGATION  WATER  TO  BE  APPLIED 
ENTER  1  FOR  YES  AND  0  FOR  NO 


ENTER  THE  FOLLOWING  RAIN  WATER  ANALYSIS 
(UNITS  ARE  MEQ/L) 

ENTER  NH4+ 


ENTER  N03- 


ENTER  CA++ 


ENTER  NA+ 


ENTER  MG++ 
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ENTER  HC03- 

0 

0 

ENTER  CL- 

0 

0 

ENTER  C03= 

0 

0 

ENTER  S04= 

0 

0 

110 


53419 


124 


30 


ENTER  THE  FOLLOWING  CROP  GROWTH  DATA 


ENTER  CROP  MATURITY  CLASS  FOR  CORN  (IN  DAYS) 
EX.  85=85  DAY  CORN,  110=110  DAY  CORN,  ETC. 


ENTER  PLANT  POPULATION  (PLANTS/HA) 


ENTER  UP  TO  10  CHARACTER  NAME  OF  CROP  VARIETY 
PIONEER 


THIS  PORTION  OF  THE  PROGRAM  GUIDES  THE  USER  IN 
ENTERING  TILLAGE  INFORMATION 

ENTER  NUMBER  OF  TILLAGE  EVENTS  (MAX=6) 


ENTER  JULIAN  DATE  OF  TILLAGE  EVENT  1 


ENTER  DEPTH  (CM)  OF  TILLAGE  FOR  EVENT  1 


ENTER  TILLAGE  CODE.   1  =  MOLDBOARD ,  1 
2  =  CHISEL  PLOW,   3  =  DISK 
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136 


0.0 


250 


323 


THE  FOLLOWING  GUIDES  IN  FERTILIZER  DATA  INPUT 
INPUT  THE  NUMBER  OF  FERTILIZER  APPLICATIONS 
ENTER  APPLICATION  DAY  NUMBER  1  (JULIAN  DATE) 


THE  FOLLOWING  WILL  SET  UP  FERTILIZER  APPLICATIONS 
FOR  EACH  FERTILIZER  APPLICATION  DATE 


ENTER  DEPTH  (CM)  OF  A  UNIFORM  FERTILIZER  APPLICATION 
FOR  APPLICATION  DATE  NUMBER     1 
(IF  =  0.0  A  SURFACE  APPLICATION  IS  INDICATED) 
DEPTH  = 


ENTER  LBS/ACRE  OF  NH4+  ADDED  IN  APPLICATION 
TYPICAL  RANGE  0.-400.  LBS/ACRE 

ENTER  LBS/ACRE  OF  N03-  ADDED  IN  APPLICATION 
TYPICAL  RANGE  0.-1600.  LBS/ACRE 

ENTER  LBS/ACRE  OF  UREA  ADDED  IN  APPLICATION 
TYPICAL  RANGE  0.-800.  LBS/ACRE 

ENTER  LBS/ACRE  OF  CA++  ADDED  IN  APPLICATION 
RANGE  0.-20000.  LBS/ACRE 

ENTER  LBS/ACRE  OF  S04=  ADDED  IN  APPLICATION 
RANGE  0.-20000.  LBS/ACRE 

ENTER  LBS/ACRE  OF  C03=  ADDED  IN  APPLICATION 
RANGE  0.-20000.  LBS/ACRE 


THE  FOLLOWING  GUIDES  IN  ORGANIC  RESIDUE  DATA  INPUT 


INPUT  NUMBER  OF  ORGANIC  APPLICATIONS 


ENTER  ORGANIC  RESIDUE  APPLICATION  DAY  (JULIAN) 
FOR  APPLICATION  NUMBER     1 


ENTER  INCORPORATION  DEPTH  (CM)  OF  UNIFORM  ORGANIC 
RESIDUE 
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0.0 


0.0 


40 


40 


0.5 


0.0 


APPLICATION     1   (MUST  BE  A  POSITIVE  NUMBER)   ENTER 
FOR  SURFACE  APPLICATION 


ENTER  CARBON: NITROGEN  RATIO  OF  ORGANIC  MATTER 
ADDED  IN  APPLICATION     1   (TYPICAL  RANGE  15-80) 
C:N  RATIO  OF  40:1  IS  ENTERED  AS  40.0 


ENTER  ORGANIC  NITROGEN  ADDED  (LBS/ACRE)  AS  OVEN 
DRY  WEIGHT   (TYPICAL  RANGE  0.-100.  LBS/ACRE) 


ENTER  DEPTH  (CM)  OF  RESIDUE  ON  SURFACE 
(TYPICAL  RANGE  0-10  CM) 


ENTER  FRESH  RESIDUE  REFLECTION  COEFFICIENT 


ENTER  RESIDUE  IDENTIFICATION  CODE 

0  IF  USER  WANTS  TO  ENTER  S,  THE  FRACTION  OF 
TOTAL  RESIDUE  WHICH  IS  LABILE 

1  IF  CORN  STOVER 

2  IF  WHEAT  STRAW 

3  IF  CELLULOSE 

4  IF  SOYBEANS 

5  IF  S  IS  TO  BE  CALCULATED  FROM  THE  EQUATION 
S  -  0.07  +  1.11(1/ACN1) ' ' (1/3) 


9  b 


Contents    of   Example 

ROSEMOUNT    1981    --    TILL   CASE 

Common   Input   File 

for   the   NTRM  Central 

.101000.    5.00.0001 

Model 

121      160           4           4           0 

0                                  1 

1        10           1           1           0 

0           0           0           0 

2           10           10 

0           111 

15.0   45.0105.0150.0 

2           14           3 

2           2           2           2 

0.0      0.0      0.0      0.0      0.0 

0.000    .1000.0000.0000.0000. 

0.10.0020.001 

0.10.0020.001 

0.10.0020.001 

0.10.0020.001 

0.10.0020.001 

0.10.0020.001 

0.10.0020.001 

0.10.0020.001 

0.10.0020.001 

0.10.0020.001 

150 


0    10 


10 


0.0   0.0   0.0   0.0 
0000.0000.0000.000 


0013.400  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00 

00   .95  0.00  0.00 

340     0.000    0.    0.   10.  500. 

0102.400  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00 


.00 

.340 

.010 

.00 

.340 

.010 

.00 

.350 

1 

1 

1 

0. 


.95  0.00  0.00 

0.000    0.    0.    5.  196. 
.600  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00 
.95  0.00  0.00 

0.000    0.    0.    1.   24. 
.600  0.00  0.00  0.00  0.00 
.95  0.00  0.00 


0.00  0.00  0.00  0.00  0.00  0.00 


0.000 
136 
323 

124  30.02 
250.    0. 
0.0  40.0  40.0 
11053419 
1.00 
200.0 
41.000   300 


0 


0. 

.5 

.50 

000 

0 


0 


0. 
0.00 


0.    0. 
1 . 0000 . 00000 . 0000 . 0000 . 000 


1       1 

0    0 
20.00 
20.00   60.00 
20.00 


20. 
11.30 

.22 
11.30 


PIONEER 
0.00 


.220 


15.30 


0.00 


1.20 


30 


7.5 


15.30 


30 


30 


08 


08 


22   20.00 


22   20.00 
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20.00      60.00 

.22 

20.00 

11.30 

15.30 

.30 

.08 

.22 

20.00 

20.00      60.00 

.22 

20.00 

11.30 

15.30 

.30 

.08 

.22 

20.00 

20.00      60.00 

.22 

Sample  Interactive     ENTER  FILE  NAME  FOR  NEW  PAN  EVAPORATION  FILE: 

Session  To  Create  Pan   NEWPAN 

Evaporation, Wind  Run,  THIS  PROGRAM  GUIDES  THE  USER  IN  THE  INPUT 

and  Solar  Radiation  OF  PAN  EVAPORATION  DATA  EITHER  ON  A  SEMI-MONTHLY  BASIS 

Input  File  OR  A  DAILY  BASIS.   IN  THE  LATTER  CASE,  DAILY  WIND 

TRAVEL  (KM/DAY)  AND  SOLAR  RADIATION  (LANGLEYS)  ARE 

ALSO  REQUIRED. 

CHOOSE  ONE  OF  THE  FOLLOWING  OPTIONS  FOR  PAN  EVAPORATION 

0.  SEMI-MONTHLY  PAN  EVAPORATION 

1.  DAILY  PAN  EVAPORATION  (+  WIND  +  SOLAR) 

1 

THE  PROGRAM  WILL  PROMPT  YOU  IN  THE  INPUT 

OF  DAILY  PAN  EVAPORATION,  WIND  TRAVEL,  AND 

SOLAR  RADIATION.   THE  PAN  EVAPORATION  IS  IN  UNITS 

OF  INCHES,  WIND  IS  IN  UNITS  OF  KM/DAY, 

AND  SOLAR  RADIATION  IS  IN  UNITS  OF  LANGLEYS/DAY. 

(EX.,  .217,193,464) 

ENTER  THE  NUMBER  OF  DAYS  OF  THE  SIMULATION  RUN  ON  THE 
NTRM  MODEL 


40 


ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY 


.10,100.0,350.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY 

.10,100.0,200.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY 

.10,100.0,300.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY 

.10,100.0,250.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY 

.10,100.0,250.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY 

.10,100.0,250.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY 

.10,100.0,250.0 
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ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY    8 

.10,100.0,250.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY    9 

.10,100.0,250.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   10 

.10,100.0,400.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   11 

.16,100.0,400.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   12 

.14,100.0,400.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   13 

.28,100.0,500.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   14 

.28,100.0,500.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   15 

.23,100.0,500.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   16 

.23,100.0,500.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   17 

.23,100.0,500.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   18 

.33,100.0,500.0 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   19 

.25,100.0,695.6 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   20 

.28,100.0,694.9 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   21 

.49,100.0,655.9 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   22 

.16,100.0,367.5 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   23 

.16,100.0,241.1 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   24 
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.16,100.0,162.2 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   25 

.12,100.0,125.4 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   26 

.08,100.0,188.1 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   27 

.10,100.0,424.4 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   28 

.10,100.0,221.8 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   29 

.26,100.0,572.8 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   30 

.26,42.5,691.3 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   31 

.26,42.5,461.3 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   32 

.39,42.5,429.9 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   33 

0.00,42.5,359.3 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   34 

.19,42.5,448.8 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   35 


.26,42.5,666.7 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   36 

.33,194.7,626.9 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   37 

.33,178.1,737.5 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   38 

.33,208.5,378.5 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   39 

0.00,218.0,298.5 

ENTER  PAN,  WIND,  AND  SOLAR  DATA  FOR  DAY   40 

.12,198.7,312.5 


99 


Contents  of  Example 
Pan,  Wind,  and  Solar 
Input  File 


.10100.0350.0 
.10100.0250.0   . 

.10100.0250.0 
.14100.0400.0 

.23100.0500.0 
.25100.0695.6 

.16100.0367.5 
.08100.0188.1 

.26100.0572.8 
0.00  42.5359.3 

.33194.7626.9 
.12198.7312.5 


.10100.0200.0 
10100.0250.0 

.10100.0250.0 
28100.0500.0 

.23100.0500.0 
28100.0694.9 

.16100.0241.1 
10100.0424.4 

.26  42.5691.3 
.19  42.5448.8 

.33178.1737.5 


.10100.0300.0  .10100.0250.0 
10100.0250.0 

.10100.0400.0  .16100.0400.0 
28100.0500.0 

.23100.0500.0  .33100.0500.0 
49100.0655.9 

.16100.0162.2  .12100.0125.4 
10100.0221.8 

.26  42.5461.3  .39  42.5429.9 
.26  42.5666.7 

.33208.5378.5  0.00218.0298.5 


Sample  Interactive 
Session  To  Create 
Precipitation  and 
Irrigation  Input  File 


ENTER  FILE  NAME  FOR  NEW  WATER  APPLICATION  FILE: 
WSAVE 


THIS  PROGRAM  GUIDES  THE  USER  IN  THE 
INPUT  OF  WATER  APPLICATIONS  DATA  FOR  NTRM 


12345 


1981 


PRECIP 


05031 


88 


ENTER  A  NUMERIC  IDENTIFIER  CODE 


ENTER  YEAR  OF  THE  RUN,  EX. ,  1981 


ENTER  NUMBER  OF  WATER  APPLICATIONS 
(MAXIMUM  NUMBER  IS  60) 


ENTER  ALPHANUMERIC  DESCRIPTION  TO  IDENTIFY 
THE  WATER  APPLICATIONS 


ENTER  MONTH, DATE  OF  MONTH, 

AND  AN  ALPHANUMERIC  IDENTIFIER  CODE 

FOR  THE  TIME  OF  WATER  APPLICATION   1  AND  TYPE 

(EX.,  0121R=JANUARY  21,  RAIN) 


ENTER  AMOUNT  OF  APPLICATION    1  AS  A 
VOLUME  PER  UNIT  AREA  (INCHES) 


ENTER  MONTH, DATE  OF  MONTH, 

AND  AN  ALPHANUMERIC  IDENTIFIER  CODE 

FOR  THE  TIME  OF  WATER  APPLICATION   2  AND  TYPE 

(EX.,  0121R=JANUARY  21,  RAIN) 


0504R 


100 


15 


0508R 


.15 


0523R 


80 


0528R 


.37 


0601R 


05 


ENTER  AMOUNT  OF  APPLICATION   2  AS  A 
VOLUME  PER  UNIT  AREA  (INCHES) 


ENTER  MONTH, DATE  OF  MONTH, 

AND  AN  ALPHANUMERIC  IDENTIFIER  CODE 

FOR  THE  TIME  OF  WATER  APPLICATION   3  AND  TYPE 

(EX.,  0121R=JANUARY  21,  RAIN) 


ENTER  AMOUNT  OF  APPLICATION    3  AS  A 
VOLUME  PER  UNIT  AREA  (INCHES) 


ENTER  MONTH, DATE  OF  MONTH, 

AND  AN  ALPHANUMERIC  IDENTIFIER  CODE 

FOR  THE  TIME  OF  WATER  APPLICATION   4  AND  TYPE 

(EX.,  0121R=JANUARY  21,  RAIN) 


ENTER  AMOUNT  OF  APPLICATION    4  AS  A 
VOLUME  PER  UNIT  AREA  (INCHES) 


ENTER  MONTH, DATE  OF  MONTH, 

AND  AN  ALPHANUMERIC  IDENTIFIER  CODE 

FOR  THE  TIME  OF  WATER  APPLICATION   5  AND  TYPE 

(EX.,  0121R=JANUARY  21,  RAIN) 


ENTER  AMOUNT  OF  APPLICATION    5  AS  A 
VOLUME  PER  UNIT  AREA  (INCHES) 


ENTER  MONTH, DATE  OF  MONTH, 

AND  AN  ALPHANUMERIC  IDENTIFIER  CODE 

FOR  THE  TIME  OF  WATER  APPLICATION   6  AND  TYPE 

(EX.,  0121R=JANUARY  21,  RAIN) 


ENTER  AMOUNT  OF  APPLICATION    6  AS  A 
VOLUME  PER  UNIT  AREA  (INCHES) 
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0602R 


.08 


0607R 


17 


ENTER  MONTH, DATE  OF  MONTH, 

AND  AN  ALPHANUMERIC  IDENTIFIER  CODE 

FOR  THE  TIME  OF  WATER  APPLICATION   7  AND  TYPE 

(EX.,  0121R=JANUARY  21,  RAIN) 


ENTER  AMOUNT  OF  APPLICATION    7  AS  A 
VOLUME  PER  UNIT  AREA  (INCHES) 


ENTER  MONTH, DATE  OF  MONTH, 

AND  AN  ALPHANUMERIC  IDENTIFIER  CODE 

FOR  THE  TIME  OF  WATER  APPLICATION   8  AND  TYPE 

(EX.,  0121R=JANUARY  21,  RAIN) 


ENTER  AMOUNT  OF  APPLICATION    8  AS  A 
VOLUME  PER  UNIT  AREA  (INCHES) 


Contents  of  Sample    12345  1981   8PRECIP 

Precipitation  and 

Irrigation  Input  05031   .88     0504R   .15     0508R   .15     0523R 

File  .80 

0528R   .37     0601R   .05     0602R   .08     0607R 
.17 
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